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CHAPTER  1 
INTRODUCTION 
1.1  Overview 

Modern  bridge  designs  often  Include  decks  whose  spans  are  at  a  skew 
angle  to  their  supports  for  economic  or  aesthetic  considerations.  The 
behavior  of  these  structures  is  largely  dependent  upon  a  variety  of 
factors  such  as  the  angle  of  skew,  direction  of  reinforcement,  aspect 
ratio,  orthotropic  stiffeners  and  types  of  loading.  The  effects  of  these 
and  other  factors  on  the  mechanics  of  deformation  may  be  interdependent, 
thus  requiring  a  study  including  variation  of  parameters.  This  report 
will  examine  the  effects  of  skew  angle  on  major  principal  moments  of  free 
spans  under  simulated  truck  loading.  Orthogonal  and  skewed  reinforcement 
will  also  be  discussed. 

The  finite  element  approach  has  been  selected  as  the  method  of 
analysis.  A  mesh  size  study  Is  undertaken  in  Chapter  2  to  examine  the 
effects  of  singularities  caused  by  concentrated  loads.  6TSTRUDL  [5] 
software  Is  used  for  the  stiffness  analysis  on  a  VAX  11-750  digital 
computer.  A  service  load  simulator  program  Is  written  in  FORTRAN-77 
(Appendix  A)  and  Is  run  semi -Interactively  with  6TSTRUDL  to  incrementally 
adjust  the  truck  position.  The  parameter  studies  are  discussed  in 
Chapter  3.  Here  the  major  principal  moment  at  center  span  is  plotted 
against  the  location  of  the  truck  for  various  angles  of  deck  skew.  The 
moment  reductions  which  are  found  for  increased  angles  of  skew  do  not 
necessarily  allow  for  a  reduction  in  reinforcing  steel.  The  angles  at 
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which  the  major  principal  moments  Intersect  the  steel  plays  an  Important 
role  In  the  efficiency  of  resisting  flexure.  This  relationship  Is 
discussed  In  Chapter  4. 

1.2  The  Nature  of  Plate  Bending 

An  orthogonal  Isotropic  plate  subjected  to  transverse  loading 
deflects  according  to  the  differential  equation  [183 

+  2  -ii L_.  +  -  -  q/D  (1) 

ix*  3x*3y*  iy* 


In  the  case  of  a  simple  span  such  as  a  bridge  deck  of  Figure  1.1(a)  the 
solution  to  (1)  must  also  satisfy  the  boundary  conditions 
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Along  Free  Edges 


where:  w  ■  displacement 

x,y  *  cartesian  coordinates 

v  «  poisson's  ratio 

q  ■  loading  function 

D  ■  Isotropic  stiffness 


(4) 

(5) 


These  equations  can  be  solved  for  many  loading  conditions  by  the  use  of 
Fourier  series.  For  skewed  plate  analysis,  the  Introduction  of  an 
oblique  coordinate  system,  such  as  that  shown  in  Figure  1.1(b)  is 


required  for  the  analysis.  After  transformation  of  the  Laplaclan, 
equation  (1)  becomes  [17] 


— (4r  +  2(1  +  2  sin2*)  +  4sin*  «  p(7,y) 

cos  ♦  3x  3x  3y  3lr3y  3x3y  3jr 

(6) 


where:  x  *  x  +  y  tan* 
y  »  y  sec* 

*  *  deck  angle  of  skew 
p(3T,  jX  =  transformed  loading  function 

Solutions  to  (6)  are  difficult  even  for  simple  loading  cases  and  are 
therefore  not  well  suited  to  extensive  parameter  studies,  although  Kennedy 
[10]  has  been  successful  with  series  solutions  on  a  computer.  Analytical 
methods  have  also  been  applied  by  Krettner  [12]  and  Lardy  [13],  Energy 
methods  have  been  used  by  Guzman  and  Lulsonl  [6]. 

1.3  Numerical  Methods 

Several  researchers  have  been  successful  in  analysis  of  skew  slabs  by 
the  finite  difference  [8],  [9]  and  finite  element  methods  [1],  [IS].  The 
growing  popularity  of  numerical  methods  for  skew  slab  analysis  may  be 
attributed  to  both  the  inadequacy  of  analytical  solutions  and  to  the 
advances  In  today's  computer  technology.  The  latter,  particularly,  has 
allowed  for  extensive  use  of  numerical  methods  in  applied  mechanics. 
Several  package  programs  such  as  STRUOL  [14]  are  available  to  the  public 
and  have  become  a  widely  accepted  method  for  analysis,  design  and  research 


The  Finite  Element  Method  was  selected  for  this  study  using  GTSTRUDL 
[5]  software  on  a  VAX  11-750  computer.  The  finite  element  approach  breaks 
up  or  "discretizes"  the  structure  to  be  analyzed  into  a  network  of  consti¬ 
tuent  elements.  In  plate  bending  each  of  these  elements  is  usually 
allowed  three  types  of  displacements  or  "degrees  of  freedom"  at  each 
corner  node  (two  orthogonal  rotations  and  one  transverse  displacement).  A 
mathematical  function  is  assumed  to  describe  the  displacement  variation 
between  the  nodes.  Then,  the  stiffness  properties  of  the  individual 
elements  can  be  developed  in  matrix  form.  If  compatability  between  the 
adjacent  elements  is  satisfied  then  the  solution  obtained  should 
"converge"  as  the  number  of  elements  is  increased.  Here  compatibility 
refers  to  the  fact  that  the  pieces  must  fit  together  and  that  all 
adjoining  elements  at  similar  nodes  must  have  corresponding  degrees  of 
freedom...  Strictly,  compatibility  is  not  completely  satisfied  for  the 
Bending  Plate  Parallelogram  (BPP)  element  which  is  used  in  this  study. 
Zienkiewicz  [20]  shows  it  is  not  possible  for  a  simple  polynomial  expres- 
slon  to  ensure  full  compatibility  when  only  one  displacement  and  two 
rotations  are  prescribed  at  the  nodes.  However,  experience  [20]  with  the 
BPP  element  shows  that  it  "converges"  to  a  good  engineering  approximation 
in  most  practical  cases. 

Although  many  different  types  of  elements  may  be  used  in  discretizing 
a  structure,  the  procedure  Is  fundamentally  the  same.  The  material 
properties  and  boundary  conditions  are  first  defined  for  the  problem.  The 
element  mesh  is  then  selected  and  the  structure  stiffness  matrix  is  formed 
from  the  known  element  stiffness  matrix.  Matrix  algegra  can  then  be  used 
to  solve  the  equation  relating  displacement,  moment,  strain  and  stress. 
Several  texts  (such  as  [3]  and  [20])  are  available  which  detail  the  finite 
element  process  and  give  many  diverse  applications. 


CHAPTER  2 


MESH  SIZE  STUDY 

2.1  Introduction 

The  first  step  in  analyzing  a  structure  by  the  finite  element 
method  Is  selecting  a  mesh  size  and  pattern.  The  mesh  is  a  very 
important  feature  of  the  study  and  must  be  selected  with  care.  In  this 
study  the  following  aspects  of  mesh  size,  type  and  pattern  were 
considered: 

1.  Accuracy  of  results 

2.  Economy  of  computer  time 

3f-  Ease  of  comparison  of  results  at  varied  angles  of  skew 

4.  Assurance  of  conforming  to  the  criterion  of  non-distorted 
elements. 

GTSTRUDL  provides  for  a  wide  range  of  elements  suitable  for  this 
study.  The  element  type  chosen  is  the  Bending  Plate  Parallelogram  (BPP) 
which  is  very  well  suited  to  skew  plate  analysis.  The  BPP  uses  a  fourth 
order  transverse  displacement  expansion  and  uses  three  degrees  of 
freedom  (one  displacement  and  two  rotations)  at  each  corner  node.  At  90 
degrees  the  BPP  Is  equivalent  to  the  more  familiar  Bending  Plate 
Rectangle  (BPR).  The  element  Is  not  considered  distorted  unless  it  Is 
skewed  to  an  interior  acute  angle  of  less  than  30  degrees. 

Consequently,  it  can  be  used  for  deck  skew  angles  of  up  to  60  degrees. 
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This  study  will  model  the  deck  as  a  free  span  between  simple 

supports.  The  thickness  Is  held  constant  at  18  inches  and  the  moments  X 

and  Y  are  released  at  the  joints  along  the  supported  edges.  Before 

proceeding  further  It  must  be  pointed  out  that  there  are  two  w^ys  to 

geometrically  "skew"  a  plate  from  a  rectangle  to  a  parallelogram,  both 

of  which  are  used  extensively  in  the  literature.  The  first  way  is  to 

keep  all  edges  the  same  length  and  allow  the  supports  to  move  closer 

together  as  the  angle  of  skew  Is  increased.  This  is  illustrated  in 

Figure  2.1(b).  The  other  widely  used  convention  is  to  maintain  a 

constant  distance  between  the  supports  and  allowing  the  free  edge  to 

"stretch"  as  a  skew  angle  is  increased.  Figure  2.1(c)  illustrates  this 

convention.  This  study  will  adopt  the  first  convention  for  the  simple 

reason  that  a  thin  strip  taken  parallel  to  the  inclined  side  will 

approximate  the  span  of  an  equivalent  simple  beam,  the  length  of  which 

will  be  held  constant  as  the  deck  skew  angle  Is  varied.  Moreover,  the 

reinforcement  Is  generally  laid  parallel  to  the  deck  edges  as  well  and 

* 

holding  these  lengths  constant  may  be  a  more  realistic  approach.  Thus 
for  all  the  examples  in  which  the  results  are  compared  between  the 
skewed  deck  and  an  orthogonal  deck,  the  geometries  shown  in  Figure 
2.1(a)  and  2.1(b)  are  Implied. 

2.2  Determination  of  Element  Size 

Before  skew  angle  can  be  investigated,  the  proper  element  size  must 
be  determined  by  a  convergence  study  of  the  solution  of  rectangular  deck 
of  dimensions  similar  to  those  used  in  practice.  Hereafter,  the  term 
"span"  will  refer  to  the  inclined  distance  (skewed  dimension)  between 


supports.  Also  the  ratio  of  width  to  span  will  be  abbreviated  HSR  where 


the  width  Is  the  distance  along  the  supported  edge.  The  WSR  parameter 
will  be  varied  throughout  the  study  In  much  the  same  way  as  the  aspect 
ratio  is  In  orthogonal  plate  analysis.  The  span  length  chosen  for  the 
mesh  size  study  was  30  feet  and  the  supported  edge  length  45  feet 
(WSR  =  1.50).  Three  mesh  sizes  were  Investigated  based  on  these 
dimensions. 

The  first  mesh  was  called  "RIGHT"  (Figure  2.2)  and  consisted  of  24 
elements  (4  In  the  span  direction  and  6  along  the  supports)  each  90 
Inches  square.  The  Idea  of  the  study  Is  to  start  with  this  coarse  mesh 
and  compare  the  results  to  increasingly  more  refined  meshes.  However, 
since  the  comparisons  are  made  based  upon  a  10  kip  concentrated  central 
load,  the  finer  meshes  will  not  be  approaching  a  finite  value  (because, 
as  Is  well  known  from  thin  plate  theory,  a  point  load  produces  Infinite 
stress)*-;  Therefore  it  is  anticipated  that  the  true  solution  will  lie 
somewhere  between  the  first  course  mesh  run  and  one  of  the  finer 
modifications.  A  more  accurate  finite  element  analysis  can  be  made 
using  the  load  area  equal  to  the  tire  contact  (imprint)  area,  and 
comparisons  can  be  made  with  the  "point  load"  mesh  runs  to  determine 
which  will  give  the  best  approximation.  The  second  mesh  was  titled 
“RTFINE"  (Figure  2.3)  and  was  given  8  span  elements  each  45"  and  10 
transverse  elements  each  54".  For  the  third  run  these  80  elements  were 
bisected  bilaterally  to  give  a  320  element  mesh  called  "RTEXFINE" 
(Figure  2.4).  As  expected  the  values  for  the  central  span  moment  under 
the  load  Increased  In  the  two  successive  finer  meshes  (see  Figure 
2.5).  The  values  went  from  3.25  k-ln/ln  (RIGHT)  to  3.86  k-in/ln 
(RTFINE)  to  4.48  k-ln/ln  (RTEXFINE). 
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TRANSVERSE  MAOOR  PRINCIPAL  MOMENT  ALONG  MIDSPAN 


2.2.1  Convergence 


The  problem  now  Is  to  determine  which  (If  any)  of  these  results  are 
close  to  the  moment  found  under  a  more  realistic  tire  load.  The  double 
wheel  contact  area  of  a  water  tanker  truck  was  studied  in  Reference  4. 
The  dimensions  are  approximated  here  to  a  rectangle  of  10"  x  24”. 
Interior  element  loads  (l.e.,  concentrated  loads  within  the  boundaries 
of  an  element)  can  be  approximated  by  placing  the  force  upon  Imaginary 
stringers  parallel  to  the  element  edges  and  using  the  resultant  reac¬ 
tions  at  the  nodal  points  as  a  new  collection  of  equivalent  loads.  This 
technique  will  be  used  later  in  the  program  SKEW  LOADER  as  a  method  for 
finding  statically  equivalent  nodal  loads  for  any  truck  position  on  a 
skewed  deck.  This  technique,  however.  Is  not  applicable  at  this  point 
as  the  objective  here  Is  to  avoid  the  singularities  caused  by  concen¬ 
trated  Juads  altogther.  Therefore  a  more  exact  approach  Is  required 
before  we  can  obtain  a  solid  basis  for  comparison  in  choosing  the 
appropriate  mesh  size. 

Since  GTSTRUDL  will  accept  uniformly  distributed  loads  only  over  an 
entire  element  (l.e.,  no  partial  element  loading),  an  extremely  fine 
localized  mesh  has  to  be  developed  to  surround  this  tire  load.  Extra 
care  must  be  taken  here  to  ensure  that  compatibility  conditions  are 
satisfied  between  element  and  that  no  elements  are  distorted. 

Distortion  Is  defined  In  GTSTRUDL  as  having  aspect  ratios  greater  than 
2.0  or  acute  angles  less  than  30  degrees.  Rectangular  elements  cannot 
be  used  In  this  localized  mesh  as  compatibility  cannot  be  maintained 
without  highly  distorting  element  aspect  ratios.  Therefore,  a 
triangular  pattern  was  carefully  assembled  for  the  localized  mesh 
(Figure  2.6)  using  GTSTRUDL  element  Bending  Plate  Hybrid  Triangle 
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(6PHT).  Due  to  the  hi yh  degree  of  correlation  between  the  data  obtained 
from  three  different  meshes  (RIGHT,  RTFINE  and  RTEXFINE)  for  areas  not 
in  the  vicinity  of  the  load,  the  RTFINE  mesh  was  chosen  to  house  the 
localized  triangular  mesh.  In  other  words.  Figure  2.6  was  Inserted  into 
the  hatched  portion  of  Figure  2.7.  The  central  portion  of  the  localized 
triangular  mesh  contains  four  rectangles  each  5“  x  12"  which  are  loaded 
with  a  uniform  pressure  of  0.04167  ksl.  This  modified  system  Is 
therefore  statically  equivalent  to  a  central  10  kip  concentrated  load  on 
the  plate.  This  mesh,  called  UTRTF1NE,  gave  a  central  span  moment  of 
4.107  k-ln/ln  which  falls  slightly  above  the  moment  given  by  the  RTFINE 
mesh  (see  Figure  2.8).  Since  the  moment  of  the  RTFINE  mesh  with  a 
concentrated  central  10  kip  load  is  only  6.2%  lower  than  this  "exact" 
value.  It  Is  considered  that  the  RTFINE  mesh  will  be  accurate  enough  for 
the  purpose  of  Investigating  the  variation  in  moment  with  deck  skew 
angle.  The  cost  in  computer  time  for  using  the  RTEXFINE  mesh  would  be 
Increased  exponentially  as  the  number  of  elements  are  quadrupled,  and 
would  as  equally  overestimate  the  moment  as  the  RTFINE  mesh  would  under¬ 
estimate  It. 


2.3  Effect  of  Aspect  Ratio 

Before  the  RTFINE  mesh  could  be  selected  as  the  proper  mesh  size 
for  the  study,  one  more  check  had  to  be  made.  If  the  results  showed 
substantial  changes  In  accuracy  when  the  aspect  ratio  of  the  plate  was 
varied,  then  this  would  Indicate  that  the  mesh  size  used  in  the  study 
would  have  to  vary  as  the  WSR  Is  changed.  A  similar  run  was  therefore 
made  to  compare  the  output  of  the  RTFINE  mesh  with  a  more  exact  value 
where  the  aspect  ratio  has  been  reduced  from  1.50  to  0.60.  The 
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respective  meshes  were  renamed  RTFINE060  and  UTWSR060  (Figure  2.9).  The 
mesh  of  Figure  2.6  was  renumbered  to  be  housed  in  the  mesh  of  Figure  2.9 
as  shown  and  comparisons  were  again  made  between  the  central  span 
moments  given  by  the  two  outputs.  The  results  showed  no  substantial 
change  in  behavior  at  the  new  aspect  ratio.  The  UTWSR060  (more  exact) 
gave  a  central  span  moment  of  5.69  k-in/in  while  the  RTFINE060  gave  5.44 
k-ln/ln,  or  4.3%  less  (see  Figure  2.10).  The  RTFINE  mesh  is  therefore 
considered  a  satisfactory  approximation  to  the  true  deck  behavior  and 
will  be  used  as  the  basis  for  investigating  the  effect  of  skew. 
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FIGURE  2.9  MESH  (a)  RTFINE060  (WITHOUT  INSERT) 
(b)  UTWSR060  (WITH  INSERT) 


FIGURE  2.10  COMPARISON  OF  RTFINE  AND  UTWSR060 


CHAPTER  3 


PARAMETER  STUDIES 
3.1  Introduction 

Now  that  the  mesh  RTFINE  has  been  selected  as  a  guide  for  the 
element  size,  the  effect  of  skew  angle  can  be  studied.  In  order  to 
obtain  a  realistic  evaluation  of  the  variation  in  maximum  moments,  the 
loading  used  must  correspond  to  the  dimensions  and  axle  weights  of  an 
actual  truck.  Furthermore,  a  consistent  convention  must  be  used  to 
establish  the  truck  position  for  varying  angles  of  deck  skew.  The 
principal  moment  at  key  points  of  the  slab  will  be  evaluated  by  GTSTRUDL 
finite  element  analysis  for  each  position  of  the  truck  as  it  moves 
incrementally  across  the  span.  Then  the  corresponding  results  for 
different  angles  of  skew  can  be  compared  when  plotted  on  the  same 
graph.  Since  the  WSR  may  affect  the  moment  variation  as  the  skew  angle 
is  changed,  it  too  will  be  treated  as  a  parameter. 

Figure  3.1  shows  the  models  to  be  studied.  The  principal  moments 
will  be  computed  for  WSR's  of  1.50,  1.00  and  0.75  while  the  skew  angle 
is  changed  from  0  to  40  degrees.  For  the  case  of  WSR  *  1.00  (rhombus), 
skew  of  20  degrees  will  also  be  Investigated.  In  all  cases  the  truck 
will  move  along  the  inclined  centerline  as  measured  from  support  to 
support.  In  the  case  of  the  rhombus,  the  effects  of  edge  loading  will 
also  be  examined  as  the  truck  is  moved  across  the  span  inset  at  48"  from 
the  left  free  edge.  In  all  cases  the  deck  is  simply  supported  at  the 
top  and  bottom  edges  and  free  along  the  left  and  right  edges.  This  will 
be  the  conventional  orientation  of  the  deck  throughout  this  report. 
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FIGURE  3.1  PARAMETER  STUDIES 


The  truck  chosen  for  the  investigation  was  the  70  kip  FDQT-SU4 
(Figure  3.2).  The  span  dimension  of  the  deck  was  selected  for 
convenience  to  be  32  feet  in  order  that  8  span  elements  of  48"  (Figure 
3.3)  may  be  used  as  in  the  RTFINE  deck.  The  RTFINE  actually  used  8  span 
elements  of  45",  though  the  element  study  showed  the  difference  should 
be  of  negligible  order.  The  transverse  dimensions  were  varied  from  48 
to  32  to  24  feet  to  allow  the  WSR  to  change  from  1.50  to  1.00  to  0.75, 
respectively.  Transverse  element  dimension  was  chosen  as  48"  for 
convenience.  Again,  these  parallelogram  element  dimensions  remain 
constant  as  the  skew  angle  of  the  deck  is  varied  (Figure  3.4(a)).  Note 
the  reference  position  of  the  skew  angle  in  Figure  3.4(a)  as  the 
complement  is  sometimes  used  In  the  literature. 

3.2  Procedure 

For  the  span  chosen  for  the  study,  the  maximum  centerline  moment  in 
an  orthogonal  deck  would  occur  under  axle  #3  (wheels  #5  and  #6  in  Figure 
3.2).  Wheel  #5  was  therefore  chosen  as  the  reference  from  which  the 
position  of  the  truck  will  be  measured.  The  truck  will  always  be 
positioned  as  If  it  were  moving  parallel  to  the  free  edges  of  the 
deck.  Since  the  program  SKEW  LOADER  (see  Appendix  A)  allows  for  input 
In  terms  of  skewed  coordinates,  the  position  of  wheel  #5  will  be  given 
in  such  a  coordinate  system.  The  X  Input  will  be  as  measured  from  the 
left  edge  of  the  deck  and  the  Y  Input  will  be  the  span  (inclined) 
distance  from  the  base  support.  For  center  loading,  the  position  of 
wheel  #5  will  be  Input  such  that  the  center  of  axle  #3  is  over  the 
transverse  centerline  of  the  deck  (see  Figure  3.4(b)  for  an  example). 

The  truck  will  be  "moved"  along  the  span  direction  while  maintaining 
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NOTE  CENTERLINE  OF  TRUCK  COINCIDES  WITH 
CENTERLINE  OF  DECK 


this  centered  position.  In  this  way  the  advantage  of  symmetry  can  be 
maintained  on  the  right  deck  and  axi symmetric  loading  can  be  maintained 
on  the  skewed  decks. 

The  following  discussion  applies  to  the  analysis  procedure  for  load 
cases  A  and  6  (see  Figure  3.1)  where  the  WSR  is  1.50  and  the  skew  angle 
Is  0  and  40  degrees,  respectively.  The  remaining  cases  are  handled  In 
an  analogous  manner  except  for  the  edge  loading  cases  which  will  be 
discussed  separately.  The  first  position  of  the  truck  Is  located  at  the 
transverse  center  and  96"  (l.e.,  2  elements)  up  from  the  base  support. 

A  finite  element  analysis  Is  then  run  with  the  truck  In  this  position 
and  the  principal  moments  at  joint  59  (deck  center)  are  obtained.  The 
truck  Is  then  moved  up  12"  (l.e.,  one  quarter  element)  In  the  span 
direction  to  a  new  position  108"  from  the  base  support,  finite  element 
analysis  . Is  again  run  with  the  truck  in  the  new  position  and  the 
principal  moments  at  joint  59  again  noted.  This  procedure  Is  repeated 
at  Increments  of  12"  until  the  truck  is  96*  (l.e.,  2  elements)  away  from 
the  far  support.  The  largest  principal  moment  at  joint  59  is  plotted 
against  the  truck  position  In  Figure  3.5  where  the  letter  label  of  the 
graph  corresponds  to  the  respective  case  In  Figure  3.1. 

The  above  procedure  Is  repeated  for  each  of  the  remaining  cases  and 
the  results  are  plotted  for  varied  angles  of  skew  at  each  WSR  to  obtain 
Figures  3.5,  3.7  and  3.9.  In  this  way  the  effects  of  skew  can  be 
Inspected  separately  for  each  WSR.  Now  for  each  of  these  graphs  the 
transverse  variation  In  moment  must  be  studied  on  neighboring  nodes  In 
order  to  be  sure  that  the  moment  at  the  transverse  center  Is  actually  at 
or  near  the  maximum  In  the  deck.  This  Is  due  to  the  fact  that  as  deck 
skew  angle  Is  Increased,  the  element  edges  are  moved  closer  to  the 


wheels  of  the  truck  which  may  allow  for  greater  localized  moments  at 
joints  off-center.  Thus  each  of  the  Figures  3.5,  3.7  and  3.9  showing 
moment  variation  with  truck  position  Is  followed  by  a  Figure  3.6,  3.8, 
and  3.10,  respectively,  showing  the  section  moment  variation  across 
midspan  In  the  vicinity  of  the  peak  moment.  This  section  diagram 
examines  the  moment  variation  at  midspan  across  two  elements  to  the  left 
and  right  of  the  peak  and  therefore  plots  for  a  total  of  5  joints.  For 
example,  the  joint  59  moment-position  graph  of  case  A  (WSR  »  1.50, 
skew  =  0.0)  peaks  when  wheel  #5  is  positioned  192"  from  the  base 
support.  Therefore  to  Investigate  the  transverse  variation  In  moment 
for  this  position  of  the  truck.  Figure  3.6  plots  the  moments  at  joints 
41,  50,  59,  68  and  77.  This  graph  is  labeled  "A-192-RT*  giving  the  load 
case,  truck  position  and  deck  skew  angle,  respectively. 

For  the  case  of  edge  loading  a  similar  procedure  Is  followed, 

though  edge  loading  Is  considered  only  for  the  cases  with  a  WSR  of 

1.00.  The  truck  Is  again  placed  96"  from  the  base  support  but  now  48" 

% 

(l.e.,  one  element)  In  from  the  left  free  edge.  A  finite  element 
analysis  Is  run  with  the  truck  In  this  position  and  the  moment  at  joint 
5  Is  noted.  The  truck  Is  then  moved  across  the  span  In  Increments  of 
12"  with  a  finite  element  analysis  run  for  each  successive  position.  In 
each  case  the  moment  at  joint  5  Is  noted,  and  as  before  the  truck  Is 
stopped  at  96"  from  the  far  support.  The  moment  at  joint  5  (the  edge 
node  at  midspan)  Is  plotted  against  each  position  of  the  truck  and  is 
shown  In  Figure  3.11  for  skew  angles  of  0,  20  and  40  degrees.  In  a 
similar  procedure  as  for  the  center  loading  cases,  the  transverse  moment 
variation  Is  plotted  In  Figure  3.12  which  examines  the  moments  across  4 
elements  starting  from  the  left  edge.  For  this  parameter  study. 


WSR  *  1.50 
CENTER  LOAD 


VARIATION  WITH  TRUCK  POSITION 


CENTER  LOAD 


WSR  •  1.00 
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[ENT  VARIATION  WITH  TRUCK  POSITION 


JOINT  NUMBER 

FIGURE  3.8  MOMENT  VARIATION  ACROSS  MIDSPAN 


CENTER  LOAD 


MOMENT  VARIATION  WITH  TRUCK  POSITION 


POSITION  OF  WHEEL  #5 

WENT  VARIATION  WITH  TRUCK  POSITION 
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Including  both  the  center  and  edge  loading  cases,  a  total  of  170  finite 


element  runs  were  made. 


3.3  Results 


In  many  of  the  skewed  decks,  the  maximum  moment  does  not  occur  at 
the  same  joint  as  in  the  corresponding  orthogonal  deck.  However,  It  Is 
seen  that  the  increase  in  moment  from  the  corresponding  joint  never 
exceeds  3%.  Therefore  the  moment  at  the  same  joint  will  be  used  for 
comparison.  It  must  be  noted  however,  that  the  moments  generally  do  not 
peak  for  the  same  position  of  the  truck.  For  example.  Figure  3.9  (HSR  * 
0.75,  center  load)  shows  that  the  right  deck  attains  maximum  moment  at 
joint  32  when  the  truck  position  (i.e.,  wheel  #5)  is  at  192"  from  the 
base  support.  The  40  degree  deck  however  peaks  earlier  when  the  truck 
position  Is  at  144".  This  is  the  general  trend  in  all  center  load  cases 
studied  here.  The  right  deck  always  peaks  at  192“  while  the  40  degree 
deck  peaks  about  48"  (one  element)  sooner.  As  might  be  expected,  the  20 
degree  deck  peaks  about  24"  earlier. 

The  fact  that  the  moments  peak  earlier  for  Increased  angles  of  skew 
can  be  attributed  to  the  effect  of  wheel  #6  (also  a  critical  wheel) 
reaching  the  deck  center  sooner.  This  allows  the  moment  at  the  deck 
center  to  peak  earlier  and  also  accounts  for  a  portion  of  the  reduction 
in  maximum  principal  moment  from  the  right  deck  because  both  wheels  5 
and  6  are  not  lying  along  the  span  center  at  the  same  time.  The  right 
deck  peaking  at  192"  Is  to  be  expected  as  this  corresponds  to  the 
position  where  the  critical  wheels  are  at  span  center.  For  the  case  of 
the  edge  loading  however,  it  Is  noted  that  the  joint  5  moment  peaks 
occur  when  the  truck  Is  positioned  at  192"  regardless  of  the  angle  of 
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skew.  Figure  3.12  shows  that  the  principal  span  moment  decreases  from 
the  left  edge  of  the  deck  at  low  angles  of  skew,  though  It  increases 
slightly  at  joint  14  at  40  degrees.  This  would  seem  to  indicate  that 
the  localized  effects  of  the  free  edge  are  less  pronounced  at  greater 
angles  of  skew. 

Table  3.1  summarizes  the  key  data  for  the  parameter  study.  It  is 
Interesting  to  note  that  the  greater  WSR's  allow  for  the  greater 
percentage  in  principal  moment  reduction  for  Increased  angles  of  skew. 
For  example,  it  is  seen  that  the  USR  of  1.50  has  a  moment  reduction  of 
24.3%  In  moving  from  the  right  deck  to  40  degree  skew  deck.  This  can  be 
compared  to  only  a  15.6%  reduction  for  the  WSR  of  0.75.  The  rate  of 
moment  reduction  is  also  non-linear  as  evidenced  by  the  table  entries  of 
WSR  =  1.00  (center  loading).  Here  the  moment  reduction  for  the  first  20 
degrees^  (i.e.,  0  to  20)  is  only  5.4%  while  for  the  second  20  degrees 
(I.e.,  20  to  40)  Is  13.5%.  Therefore  the  rate  of  moment  reduction  Is 
increased  with  an  increase  In  skew  angle.  In  the  case  of  edge  loading, 
the  behavior  Is  similar  except  that  the  percentages  are  somewhat 
Increased.  It  Is  Important  to  note  that  these  decreases  in  maximum 
principal  moment  do  not  necessarily  allow  for  commensurate  reduction  in 
reinforcing  steel.  The  angle  at  which  the  reinforcement  Is  laid  and  the 
direction  of  these  principal  moments  play  an  Important  role  in  the 
amount  of  steel  required.  These  effects  will  be  discussed  In  greater 
detail  in  Chapter  4. 

3.4  Contour  Plot  Description 


Appendix  C  shows  the  contours  of  major  principal  moment  and  Z 
(transverse)  displacement  for  skew  angles  of  0  and  40  degrees.  These 
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plots  were  obtained  from  GTSTRUDL  graphics  and  plotted  on  a  Tektronix 
print  device  from  the  scope  environment.*  The  loading  conditions 
selected  for  comparison  were  a  central  10  kip  load  and  the  peak  SU-4 
central  load.  The  WSR  Is  1.50  throughout. 

Figures  C.l  and  C.2  show  the  effect  of  skew  angle  on  the  major 
principal  moments  under  a  central  10  kip  load.  Note  that  the  contour 
gradient  remains  essentially  perpendicular  to  the  supports.  This  same 
effect  can  be  seen  on  Figures  C.3  and  C.4  (SU-4  peak  loading  case)  where 
the  path  of  "steepest  descent"  Is  nearly  along  the  shortest  line  to  the 
supports.  The  concentration  of  contours  near  the  obtuse  corner  show 
that  the  major  principal  moment  Increases  more  quickly  (along  a  line 
towards  the  center)  than  at  the  acute  corner.  Figures  C.3  and  C.4  can 
be  compared  to  the  transverse  moment  variation  graphs  for  load  cases  A 
and  B  as  shown  In  Figure  3.6. 

Figures  C.5  through  C.8  Illustrate  the  variation  In  deflected  shape 
between  skew  angles  of  0  and  40  degrees.  Figures  C.4  and  C.5  refer  to 
the  central  10  kip  load  while  C.7  and  C.8  refer  to  the  SU-4  peak  load. 

It  Is  seen  that  the  behavior  Is  similar  for  both  loading  cases.  As 
would  be  expected  physically,  the  displacement  gradient  Is  perpendicular 
to  the  supports  and  the  direction  of  "minimum  descent”  Is  along  the  span 
centerline. 
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1  The  scope  environment  Is  a  characteristic  operating  domain  In 

GTSTRUDL  which  allows  for  Interactive  graphics  at  the  terminal  CRT. 


CHAPTER  4 


REINFORCEMENT 
4.1  Introduction 

The  safe  and  economic  proportioning  of  reinforcement  Is  critical  in 
bridge  deck  design.  Ideally,  the  re-bars  should  be  laid  orthogonal  to 
the  principal  design  moments.  However  since  the  principal  angles  vary 
from  point  to  point  throughout  the  deck  for  even  a  single  loading  case, 
the  reinforcement  cannot  be  placed  ideally  In  a  practical  sense. 
Furthermore  since  design  is  often  based  on  a  series  of  different  loading 
cases,  the  principal  angle  will  often  vary  at  the  point  as  well. 
Therefore  the  concept  of  ideal  reinforcement  for  a  bridge  deck  Is  a 
trivial  one. 

Certain  general  directions  for  the  reinforcement  are  a  good  deal 
more  economical  (in  terms  of  required  steel  quantity  only)  than 
others.  Cope  [2]  has  studied  orthogonal  reinforcement  for  skewed  decks 
and  compared  the  experimental  results  of  placement  parallel  to  the 
supports  versus  parallel  to  the  skewed  edge.  The  results  of  this  study 
will  be  discussed  in  more  detail  in  section  4.3.  Orthogonal 
reinforcement  in  skewed  decks  however,  has  severe  limitations  In 
practice.  Since  a  large  number  of  re -bars  of  different  lengths  are 
required  for  this  design,  a  great  deal  of  extra  labor  is  required  for 
cutting  and  placing  the  reinforcement.  Therefore,  though  It  Is 
generally  deslreable  to  have  reinforcing  steel  running  perpendicular  to 
the  supports  (to  resist  the  span  moment),  such  an  arrangement  is  usually 
not  practical. 


Discussions  with  FDOT  engineers  Indicate  that  most  reinforcement  in 


skewed  decks  today  is  not  orthogonal.  The  primary  reinforcement  is 
generally  laid  parallel  to  the  deck  edges  and  the  transverse  reinforce¬ 
ment  is  laid  parallel  to  the  supports.  This  is  certainly  the  easiest 
from  a  construction  point  of  view,  but  for  largely  skewed  bridges  it  may 
require  a  considerable  amount  of  extra  reinforcement.  Morely  [16]  has 
developed  a  design  process  for  non -orthogonal  reinforcement  in  skewed 
decks.  Though  the  method  may  be  tedious  for  practical  design  problems 
witout  the  aid  of  computer,  it  serves  well  to  illustrate  the  effects  of 
skewed  reinforcement.  This  design  will  be  discussed  In  section  4.4  and 
some  examples  worked  which  are  relevant  to  the  parameter  study  of 
Chapter  3. 

Even  though  the  STRUDL  analyses  are  purely  elastic.  It  Is  relevant 

to  first. examine  the  failure  mechanisms  of  skewed  decks  by  the  theory  of 

yield  lines  before  discussing  reinforcement  in  detail.  For  a  more 

complete  discussion  on  yield  line  theory  and  applications,  several  texts 

* 

are  available  such  as  [7],  [19]. 

4.2  Failure  Analysis  by  Yield  Lines 

Analysis  of  plate  capacity  can  be  determined  theoretically  by  the 
principle  of  virtual  work  and  an  assumption  of  a  yield  line  pattern. 

The  work  done  externally  by  the  applied  forces  acting  through  the  plate 
deflection  Is  equated  to  the  internal  work  done  by  the  rotation  of  the 
moments  acting  on  the  yield  lines.  The  collapse  load  can  then  be  solved 
for  in  terms  of  the  ultimate  unit  moment  capacity  of  the  plate.  Often 
the  main  problem  in  the  yield  line  approach  Is  determining  the  correct 
(or  critical)  failure  mechanism  to  be  analyzed.  Yield  line  patterns  can 


V. 


be  assumed  from  some  experience  or  from  some  generally  established 
guidelines.  Hughes  [7]  gives  some  such  guidelines,  although  generally 
the  critical  pattern  is  not  easily  found  even  for  some  rather  simple 
geometries. 

For  example,  Hughes  [7]  shows  three  possible  failure  mechanisms  for 
a  simply  supported  skewed  span  as  shown  in  Figure  4.1.  The  deck  is 
given  a  central  point  load  W  and  the  orthotropic  rei nf orcement  is 
oriented  parallel  to  the  deck  edges  as  is  generally  used  In  practice 
(see  Figure  4.1).  The  figure  and  equations  have  been  modified  to  agree 
with  the  conventions  used  in  this  report  thusfar.  It  Is  seen  that  the 
elliptical  fan  pattern  represents  the  least  allowable  failure  load. 

Since  bridge  decks  are  generally  not  subjected  to  a  single  central 
concentrated  load,  this  case  may  be  of  purely  academic  interest. 

However,  in  deriving  the  critical  load  for  the  elliptical  failure 
mechanism,  Hughes  outlines  (and  proves)  a  very  useful  principle  for 
obtaining  an  affine  isotropic  right  slab  from  an  orthotropic  skewed 
one.  This  procedure  is  discussed  below. 

A  skewed  deck  such  as  that  shown  in  Figure  4.2(a)  can  be 
transformed  into  an  equivalent  orthogonal  deck  for  the  purposes  of  . 
analysis.  The  equivalent  right  deck  (or  so  called  "affine  deck")  shown 
in  Figure  4.2(b)  is  of  course  much  easier  to  analyze  and  can  be 
obtained  by  the  following  rules: 

(a)  Deflections  are  Identical  at  corresponding  locations  in  both 
decks. 

(b)  Given  that  m  and  urn  are  the  ultimate  resisting  moments  in  the 
reinforcing  direction  of  the  actual  deck,  then  the  affine  deck 
has  ultimate  resisting  isotropic  moment  m. 
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(c)  Given  that  the  Y  coordinate  axis  Is  along  the  m  reinforcement 
direction  in  the  actual  deck,  then  all  distances  measured  In 
the  Y  direction  are  the  same  for  both  decks. 

(d)  Given  that  the  X  coordinate  axis  is  along  the  urn  reinforcement 
direction  in  the  actual  deck,  then  it  is  taken  at  a  right 
angle  to  the  Y  axis  for  the  affine  deck. 

(e)  In  order  to  obtain  a  dimension  In  the  X  direction  of  the 
affine  deck,  divide  the  corresponding  length  on  the  actual 
deck  by^yT 

(f)  Divide  the  loads  on  the  actual  deck  b yyfv  cos<j>  to  obtain  the 
corresponding  loads  on  the  affine  deck. 

Using  the  above  procedure,  an  interesting  result  is  obtained  if  the 
deck  in  Figure  4.2(a)  is  analyzed  assuming  a  yield  line  occurs  across 
midspan  (see  Figure  4.3).  This  assumption  for  the  failure  mechanism  may 
be  valid  for  the  multiple  load  cases  arising  on  bridge  spans.  The  load 
P  is  a  sum  of  the  concentrated  wheel  loads  acting  across  midspan  and 
should  cause  failure  to  occur  in  beam  action.  If  m  represents  the  unit 
moment  resisting  capacity  of  the  affine  deck  In  Figure  4.3(b),  then  by 
statics 


&  (jTTos^  9 
or 

m  =  4b"  cosT 

Therefore  the  deck  of  Figure  4.3(a)  is  equivalent  to  an  orthogonal  deck 
of  span  length  L  and  base  of  a  length  diminished  by  the  cosine  of  the 
angle  of  skew.  This  illustrated  in  Figure  4.4. 


EQUIVALENT  DECKS  IN  BEAM  ACTION  FAILURE 


4.3  Orthogonal  Reinforcement  in  Skewed  Decks 

An  experimental  investigation  was  conducted  by  Cope  [2]  on  the 
effects  of  orientation  of  orthogonal  reinforcement  in  45  degree  skewed 
decks.  Although  orthogonal  reinforcement  is  generally  not  used  in 
skewed  bridge  design,  the  study  finds  some  interesting  results  which  may 
lead  to  a  better  understanding  of  flexural  behavior.  The  test  procedure 
was  fundamentally  a  comparison  between  two  45  degree  skew  slabs,  each 
with  orthogonal  reinforcement.  Slab  A  had  secondary  reinforcement 
parallel  to  the  supports.  Slab  B  had  primary  reinforcement  parallel  to 
the  skewed  edges.  The  propagation  of  cracking  and  deflections  in  the 
two  slabs  were  examined  for  increasing  loads  in  various  locations. 

The  tests  showed  considerable  behavioral  differences  in  the  two 


slabs  in  the  areas  of  deflections,  cracking  and  failure  mechanisms. 
Sagging  cracks  were  first  Initialized  on  slab  A  running  towards  the  free 
edges  and  parallel  to  the  supports,  and  under  heavier  load  hogging 
cracks  appeared  at  the  obtuse  corners.  Initial  sagging  cracks  in  slat  B 
occured  earlier  (at  about  80*  of  the  load  for  A)  as  did  the  hogging 
cracks  at  the  obtuse  corners  (at  about  67%  of  the  load  for  A).  The 
cracking  patterns  were  similar  in  direction  and  spacing,  though  somewhat 
stralghter  and  more  continuous  In  slab  B.  The  modes  fo  failure  of  the 
two  decks  were  different  as  well.  Slab  A  Initially  failed  in  shear  on 
the  free  edge  at  the  obtuse  corner.  At  about  40%  higher  load,  slab  B 
developed  a  top  surface  crack  between  the  two  obtuse  corner  supports 
which  was  wide  enough  to  produce  a  discontinuity  In  slope.  Both  slabs 
were  able  to  carry  considerably  more  load  after  their  initial 
failures.  Ultimate  failures  were  punching  shear  at  the  obtuse  corner 
for  slab  A  and  excessive  deflections  for  slab  B. 


The  report  concluded  that  the  behavior  of  skewed  decks  Is  strongly 
influenced  by  the  direction  of  the  reinforcement.  When  the  reinforcing 
is  placed  orthogonally  and  parallel  to  the  supports,  the  slab  is  stiff 
and  behaves  well  under  service  loads.  However,  heavy  concentrated  loads 
are  not  distributed  well  across  the  slab  and  a  large  reactive  force 
develops  at  the  obtuse  corner.  This  led  to  local  failure  in  the  test. 
For  slabs  with  orthogonal  reinforcement  parallel  to  the  free  edges,  a 
more  flexible  slab  results  which  better  distributes  moments  due  to 
concentrated  loads.  Greater  deflections  can  occur  and  ultimate  failure 
load  is  increased,  though  again  hogging  cracks  at  the  obtuse  corner  may 
limit  servlclbllity. 

The  results  of  this  study  are  consistent  with  a  finding  by  Kennedy 
[11]  which  experimentally  Investigated  the  stress  near  corners  of  simply 
spanned  .-Skewed  plates.  Here  too  it  was  found  that  stresses  near  the 
obtuse  corners  are  significant  and  Increase  with  increasing  angles  of 
skew.  Furthermore,  the  stress  at  the  obtuse  corners  may  exceed  the 
maximum  stress  at  center  span.  Kennedy  then  recommends  that  the  obtuse 
corners  of  concrete  skewed  decks  be  heavily  reinforced  top  and  bottom  in 
directions  parallel  and  perpendicular  to  the  supports  as  well  as 
parallel  to  the  free  edge.  The  accute  corners  should  be  heavily 
reinforced  at  the  bottom  and  nominally  reinforced  at  the  top.  The 
direction  of  reinforcement  near  the  acute  corners  should  be 
perpendicular  to  a  line  between  them  and  parallel  to  the  free  edge. 

4.4  Design  for  Skewed  Steel 

A  design  procedure  was  worked  out  by  Morely  [16]  on  proportioning 
skewed  reinforcement  which  is  laid  parallel  to  the  edges  of  the  deck. 


Since  this  is  generally  the  orientation  used  in  construction,  this 
method  could  prove  useful  in  design  practice.  Equations  are  developed 
based  on  the  optimum  proportionment  of  steel  to  resist  a  moment  triad  at 
a  point  in  the  deck.  Therefore  the  method  begins  with  the  assumption 
that  elastic  analysis  has  been  performed  for  a  load  case  on  the  deck  and 
that  a  set  of  three  moments  (two  bending  and  one  torsional)  are  known  at 
the  point  considered.  The  design  procedure  is  then  carried  out  using  a 
series  of  charts,  tables  and  equations  given  in  the  reference. 

Morley's  design  for  skew  reinforcement  to  resist  a  single  moment 
triad  consists  of  the  following  steps: 

(a)  Orient  the  sign  convention  of  the  known  moment  triad  in 


accordance  with  Figure  4.5 


approximation  for  la.  (k*0  If  no  minimum  steel  is  specified) 

(c)  Enter  the  charts  with  the  nearest  values  of  4  and  k  and  locate 
the  appropriate  region 

(d)  Evaluate  alx,  a^,  a2x,  a2^»  8$  and  6h  from  the  appropriate 
equations  In  the  tables 

(e)  Compute  la  by  an  Iteration  equation  of  section  equilibrium 

(f)  Compute  Alx,  A1+,  A2x,  A^,  from  alx  -  «ylaAix  (x  ♦♦  ♦) 
where  these  variables  and  others  used  In  the  design  examples  are  defined 
as  follows: 

A^  specified  minimum  steel  are  per  unit  slab  width 

ajx,  a^  'area  functions'  with  dimensions  of  moment 

de,  d0  depths  for  bottom  and  top  steel,  respectively 

k  failure  factor  connecting  Am^n  and  applied 

twislng  moment  Mxy 


REINFORCEMENT 

DIRECTIONS 


1a  lever  arm  of  tensile  steel  In  failure  direction 

Mx,  My,  M^y  bending  and  twisting  moments  per  unit  slab  width 

Mj#  M2  principal  moments 

Tj,  T2  total  resolved  steel  forces  in  failure  direction 

at  bottom  and  top,  respectively 
e$,  8h  optimum  failure  directions  for  sagging  and 

hogging 

+  skew  angle  of  reinforcement 

Oy  yield  stress  of  steel 

0  value  of  concrete  stress  block  at  failure 

V 

A  total  of  five  examples  will  be  worked  here  to  examine  the  effects 
of  skewed  reinforcing  steel.  The  first  three  examples  show  the  extra 
reinforcement  required  for  a  20  degree  skew  deck  over  a  right  deck  when 
designing  for  the  central  span  moment  for  a  100  kip  load.  Examples  4 
and  5  similarly  illustrate  for  the  peak  moments  produced  by  SU-4 
truck.  It  should  be  noted  that  this  design  procedure  could  prove  very 
tedious  for  practical  problems  without  the  aid  of  a  computer.  This  Is 
because  the  moments  at  many  points  in  the  slab  would  have  to  be 
considered  for  various  combinations  and  positions  of  vehicle  loads.  The 
results  for  each  trial  run  are  very  sensitive  to  the  sign  and  magnitude 
of  the  torsional  moment  at  the  point  considered.  These  multiple  moment 
triads  must  be  considered  and  the  minuimum  steel  retained  at  each 
successive  Iteration.  Morley  makes  allowances  for  this  In  an  extension 
of  the  procedure  for  multiple  moment  triads  along  with  some  guidelines 
for  a  rational  approach  to  design. 

It  should  be  understood  that  In  the  examples  which  follow,  the 
steel  in  the  slab  is  not  completely  designed  for  the  given  loading 


condition,  since  only  the  moment  at  the  geometric  center  Is  considered  In 
each  case.  The  examples  are  Intended  only  to  Illustrate  the  effects  of 
non-orthogonal  reinforcement  at  a  point  for  a  few  different  load  condi¬ 
tions.  Morel y  emphasizes  the  effects  of  the  sign  of  the  torsional  moment 
by  working  two  similar  examples  of  a  20  degree  skew  deck.  The  first  has 
a  point  subjected  to  My  »  +50  k-ln/ln,  Mx  ■  +35  k-ln/ln  and  MXy  =  +20 
k-ln/ln  and  the  total  required  steel  Is  found  to  be  3.96  in2/ft.  In  the 
second  example  the  sign  of  MXy  Is  changed  to  -20  k-ln/ln  and  the  total 
amount  of  steel  becomes  8.20  in2/ft,  or  roughly  twice  the  amount.  This 
Increase  occurs  because  the  direction  of  the  major  principal  moment  Is 
shifted,  due  to  the  sign  change  of  MXy,  from  within  the  acute  angle 
between  the  reinforcement  to  the  obtuse  angle.  This  reduces  the  unit 
effective  relstance  of  the  steel,  and  hence  a  greater  amount  Is  required. 

For  the  examples  discussed  below,  reference  Is  made  to  Figure  4.6 
and  the  following  data  Is  applicable  to  each: 

Deck  thickness  =  18.0  In 
E  «  3324  ksl 
f'c  -  3401  psl 
o  •  2891  psl 

v 

VISR  *  1.00 
Oy  *  50  ksl 

In  addition,  the  following  equations  from  Morley  [16]  are  common  to  all 
examples: 

1.  Tj  -  y-  [alx  cos2e  +  a^  s1nZ(e  -  ♦)] 

T«  fT.  +  Tp\2 
>'•  \m\*  (T7>  do  -  2~~o~~Tj 


FIGURE  4.6  SKEWED  REINFORCEMENT  EXAMPLES 


Example  1.  Right  deck  100  kip  central  load 


Step  (a)  GTSTRUDL  gives  the  following  moments  in  the  appropriate  sign 
convention 

Mx  «  +  28.49  k-in/ln 
My  *  +  43.03  k -in/in 
MXy  *  0.0  k-in/in 


(b) 


\ 

TV 


M. 


IK 


xy 1 


(c)  Using  Figure  2  [16]  the  appropriate  region  Is  A 

(d)  From  Table  1  [16] 

a2x  *  a2y  *  0  ■*•  No  Top  Steel  Required 
alx  *  Mx  +  lMxy  I  =  28.49 
aiy  =  My  +  |MXy I  -  43.03 

®s  "  ”45° 

(e)  From  equation  I,  Ti  =  [28.49  cos2(-45)  +  43.03  sin2(-45)] 

a 

Tt  «  35.76/1 a  ,  T2  *  0 

T, 

From  equation  II,  la  «  16  - 

or 

la2  -  16  la  +  6.18  «  0  ♦  la  - 

(f)  From  equation  III,  12(28.49)  ■  50(15.60)Alx 

Alx  *  r*44  l"Z/ft 
12(43.03)  -  50(15.60)Aly 
Aly  -  0.66  1n2/ft 


15.60 


59 


!•  *. 

i 

I 

0 

L 


Example  2.  20  deg  skew  deck  using  moment  triad  of  example  1 


Step  (a) 


’xy 


+  28.49  k-in/in 
+  43.03  k-in/in 
0.0  k-in/in 


Mx  My 

(b)  ■Hy’TV*" 

(c)  Using  Figure  3  [16]  the  appropriate  region  Is  A2 

(d)  From  Table  1  [16] 

a2x  *  a2+  “  0  ^  Top  Steel  Re<lu1red 

•lx  *  "x  *  "y  t»"2*  +  [My  **"♦  •  V 

alx  »  28.49  -  43.03(tan  20)2  +  l1  t43-03 

alx  *  50.85 

*U  '  ^  +  25TT  |My  un*  ’  V 


43.03 


U  (cos  20) 
aj^  *  65.36 

°s  *  7  (90  +  ♦) 
e  -  +  55° 


7  +  cos"20  '  I43*03  tan  20  '  °l 


(e)  From  equation  I,  ■  -y —  [50.85  (cos  20)2  +  65.36 


Tx  -  66 .4/1 a,  T2  »  0 


From  equation  II,  la  *  16  -  ^-789lT 

or 

la2  -  161 a  +  11.48  -  0 
1.  «  15.25" 


t 


tan  20-0- 


(sin  35)2- 


I 


(f)  From  equation  III,  12(50.85)  *  50(15.25)  Alx 

Alx  -  0.80  in2/ft 
12(65.36)  *  50(15. 25)Alf 
A^  *  1.03  in2/ft 

Comparison  of  this  result  with  Example  1  Illustrates  the  effects  of  20  deg 
skew  reinforcement  minimized  for  the  same  moment  triad.  The  total  steel 
per  unit  width  becomes  (0.80  +  1.03)  cos  20  =  1.72  in2/ft  or  56%  more  than 
Example  1. 

Example  3.  20  deg  skew  deck  100  kip  central  load 


Step  (a)  GTSTRUDL  gives  the  following  moments  In  the  appropriate  sign 
convention 


*-jMx  =  +  24.36  k-in/in 
My  *  +  40.77  k -In/in 
MXy  *  +  0.646  k-ln/ln 


37-71  T^r 


40.77 

rare 


63.11 


(c)  Using  Figure  3  [16]  the  appropriate  region  Is  A2 

(d)  From  Table  1  [16] 

*2x  *  a24  *  0  *  N°  top  steel  required 

»ix  -  Mx  - ", ta"2*  ♦  L1-a i  [M, t,n  ♦  •  Hx,l 

alx  -  24.36  -  40.77  (tan  20)2  +  C40*77  tan  20-0. 

ajx  ■  44.39 

■l#  -^?  +  coTT 


y  v 


\K  tan  *  - 


*u  *  -  +  75 rar  I40*77  tan  20  *  °*6461 

14  (cos  20/  cos  20 
a1+  »  61.24 

6s  58  +  7  (90  +  ♦) 


e  *  +  55° 

$ 

(e)  From  equation  I,  T\  »  <J~  [44.39  (cos  55)2  +  61.24(sin  35)2] 

a 

Tj  »  34.7/la  ,  T2  -  0 

From  equation  II,  la  *  16  - 

or 

la2  *  -161 a  +  6.0  =  0 
1a  *  15.62" 

(f)  From  equation  III,  12(44.39)  =  50(15. 62) Ajx 

Alx  =  0.68  1n2/ft 
12(61.24)  *  50(15. 62) Aj^ 

Aj^  «  0.94  in2/ft 

Comparison  of  this  result  with  Example  1  Illustrates  the  additional 
steel  required  at  the  geometric  center  for  the  20  deg  skew  deck  vs.  the 
right  deck,  both  under  a  central  100  kip  load.  The  total  steel  per  unit 
width  becomes  (0.68  +  0.94)  cos  20  *  1.52  1n2/ft  or  38%  more  than 
Example  1.  This  Increase  Is  due  to  the  major  principal  moment  lying 
within  the  obtuse  angle  of  the  reinforcement  as  shown  In  Figure  4.7. 


Example  4.  Right  deck  peak  SU-4  loading 


Step  (a)  6TSTRUDL  gives  the  following  moments  in  the  appropriate  sign 
convention 

Mx  =  +  5.07  k-in/in 
My  =  +  13.90  k-in/in 
MXy  =  0.0  k-ln/ln 

H  M 

(b)  T«^*<T** 


(c)  Using  Figure  2  [163  the  appropriate  region  Is  A 

(d)  From  Table  1  [16] 

a2x  =  a2y  “  0  +  No  Top  Steel  Required 

“~alx  8  Mx  +  lMxyl  8  5*07 
aly  -  My  +  |Mxy|  -  13.90 
e  *  -45° 

(e)  From  equation  I,  Tj  *  9.49/la,  Tg  *  0 


From  equation  II,  la 

or 


16  - 


2(2.891) 


la2  *  -161 a  *  1.64  -  0 
la  =  15.90" 

(f)  From  equation  III,  12(5.07)  *  50(15.90)Ajx 

Alx  =  0.08  in2/ft 
12(13.90)  -  50(15.9Q)Aly 
Aiw  s  0.21  in2/ft 


vs 


64 


Example  5.  20  deg  skew  deck  peak  SU-4  loading 


mm 


Step  (a)  GTSTRUDL  gives  the  following  moments  in  the  appropriate  sign 
convention 

Mx  *  +4.66  k-ln/in 
My  *  +13.13  k-ln/in 
Mxy  =  -0.36  k-ln/in 


/.i  x  _  4.66  _  10  n  y  .  13.13  _  E 

(b)  TCJ  Off  12’9  TTTT  -  05“  36*5 


l  . 


(c)  Using  Figure  3  [16]  and  noting  that  Mxy  <  0 

Mx  ♦  -4.66  ,  My  -13.13 
and  the  appropriate  region  is  B 

(d)  From  Table  1  [16]  with  suffixes  1,  2,  s,  h  Interchanged 
a2x  *  a2*  =  0  +  No  T°P  Steel  Required 


"  "Mx  +  My  tan2*  -  — ToT-*1--*-  ("»  tan  ♦  -  ",  J 


alx  "  "Mx  +  My 


alx  «  4.66  -  13.13  {tan  20)2  -  ^1^Q2-|g-n--2-^  [-13.13  tan  20  ♦  0.36] 


alx  *  10.84 

M.. 


‘M  '  3*  *  |M>  *“  *  '  V 

'w  *  -~^r +  corro  I-13-13  20  *  °-36l 


au  -  19.57 


es  -  ♦  ^  (90  +  ♦)  ■  55° 


(e)  From  equation  I, 


Tx  =  1  [10.84  (cos  20) 2  +  19.57  (sin  35)2] 
a 

Ti  =  16.0/la,  T2  *  0 


From  equation  II,  la 

or 


16  - 


2(2.891) 


la2  ~  161a  +  2-77  ■  0  * 
(f)  From  equation  III,  12(10.84)  *  50(15. 82)Ajx 

Alx  «  0.16  1n2/ft 
12(19.57)  *  50(15.82)AU 
A1+  -  0.30  In2/ ft 


15.82 


Comparison  of  this  result  with  Example  4  Illustrates  the  additional 
steel  required  at  the  geometric  center  for  the  20  deg  skew  deck  vs.  the 
right  deck,  both  under  peak  SU-4  loading.  The  total  steel  per  unit 
width  ..becomes  (0.16  +  0.30)  cos  20  =  0.43  1n2/ft  or  49%  more  than 


Example  4.  This  increase  is  again  due  to  the  major  principal  moment 
lying  in  the  obtuse  corner  of  the  reinforcement  as  shown  in  Figure  4.8 


3  (k-in/in) 


CHAPTER  5 


CONCLUSIONS 


5.1  Summary 


’The  finite  element  method  was  used  to  examine 


(thPeff 


effects  of  skew 


angle  on  the  major  principal  moments  of  simple  spans.  A  mesh  size  of 
eight  span  elements  in  30  ^fe^Twas  found  to  give  results  within  about  6% 
of  those  found  from  a  more  intricate  localized  mesh  surrounding  a  more 
realistic  tire  load.  These  results  varied  little  with  change  in  aspect 
ratio.  This  mesh  size  was  therefore  chosen  as  a  guide  for  the  parameter 


studies  which  followed. 


wer 


Skev*  angles  of  0,  20  and  40  degreearwere  investigated  for  various 
width  to  span  ratios  (WSR's)t  6s  shown  in- Figure  The  major 

principal  moments  at  key  points  in  the  deck  were  examined  as  an  FDOT 


SU-4  type  truck  was  moved  across  the  span.  The  service  load  simulator 
program  tfTAppenU  I  a  A)>was  used  to  calculate  the  equivalent  nodal  forces 
for  the  truck  In  any  position  on  the  skewed  deck,  fao-^esults  showed 
that  for  center  loading,  peak  moment  reductions  of  up  to  24%  were  found 


for  40 


deg&e^ske 


skew  decks  over  orthogonal  decks 


.  ^Th^peri 


percent  reductions 


were  lower  for  smaller  WSR's.  In  general,  the  major  principal  moments 
peaked  earlier  (i.e.,  at  lesser  advanced  positions  of  the  truck  along 
the  span)  for  Increased  angles  of  skew^This  can  be  attributed  to  the 
fact  that  the  wheels  of  the  axles  do  not  reach  span  center 


simultaneously 


>The  direction  of  reinforcement  plays  an  important  role  in  the 
flexural  behavior  of  the  deck .tfclJ-^Jrthogonal  reinforcement  laid 
parallel  to  the  supports  provides  fd*  a  stiffer  slab  which  behaves  well 
under  service  loads.  Orthogonal  reinforcement  laid  parallel  to  the  free 
edges  provides  a  more  flexible  slab.  HoweverT^kewed  reinforcement  Is 
the  preferred  method  from  the  construction  point  of  view.  The  affine 
deck  method  discussed  in  section  4.2  seems  to  be  the  most  practical  for 
use  in  the  design  office.  Horely's  [16]  method  for  minimum  steel  at  a 
point  would  seem  to  have  limitations  in  practice,  although  it  serves 
well  to  Illustrate  the  considerable  effect  that  the  direction  of  the 
major  principal  moment  has  on  the  quantity  of  skewed  steel  required. 


5.2  Further  Study 


Metre,  extensive  parameter  studies  could  be  used  with  different  truck 
types  and  combination  loading.  Program  SKEW  LOADER  (Appendix  A)  could 
be  easily  modified  for  this  application.  In  addition,  the  effects  of 
girders,  deck  thickness  and  material  properties  could  also  be  studied. 
The  reinforcement  design  method  developed  by  Morely  [16]  is  well  suited 
to  computer  programming.  This  could  be  an  area  of  further  study  along 
with  complete  design  examples  and  comparison  with  the  affine  deck 
method. 

The  stress  concentration  effect  at  the  obtuse  corners  could  be 
investigated  by  the  finite  element  method.  This  could  be  of  particular 
Importance  as  experiments  have  shown  that  this  is  often  the  area  of 
Initial  failure.  An  extremely  fine  local  mesh  could  be  assembled  (such 
as  along  the  line  of  Figure  2.6)  in  the  area  of  the  obtuse  corner  and 
the  moments  Investigated  with  variation  in  deck  angle  of  skew. 


APPENDIX  A 


PROGRAM  SKEW  LOADER 
A.l  Introduction 

Since  the  parameter  study  calls  for  moving  a  truck  Incrementally 
across  a  skewed  deck  and  evaluating  the  moments  for  each  truck  position, 
a  method  had  to  be  developed  for  loading  the  deck.  CTSTRUDL  allows  for 
concentrated  loads  to  be  placed  only  at  the  nodal  points  of  the  mesh. 
Thus  the  weight  of  each  truck  wheel  within  an  element  must  be  broken 
Into  Its  statically  equivalent  nodal  forces  before  It  can  be  Input  to 
GTSTRUDL  for  stiffness  analysis.  To  do  this  by  hand  would  be  extremely 

tedious  as  4  reactions  would  have  to  be  calculated  for  each  of  8  wheels 

for  17(f  different  load  cases.  Furthermore,  the  positions  of  the  wheels 
In  relation  to  the  mesh  nodes  are  geometrically  complex  making  large 
numbers  of  calculations  prohibitive.  This  type  of  problem  Is  therefore 
well  suited  to  computer  programming. 

Program  SKEW  LOADER  was  developed  to  output  the  equivalent  nodal 
forces  for  a  truck  positioned  on  a  skewed  deck.  The  program  will 
generate  any  mesh  size  on  a  deck  of  any  dimensions  at  any  angle  of  skew. 

The  data  for  the  truck  FOOT  SU-4  Is  stored  Internally  In  the  program. 

Although  any  type  of  truck  may  be  used  If  the  truck  width,  axle  spacing 
and  axle  weights  are  known.  The  user  also  has  the  option  of  specifying 
the  truck  position  In  rectangular  or  skewed  coordinates.  The  latter  is 
recommended  for  most  applications  and  In  general  greatly  simplifies  the 
Input  procedure.  The  position  of  any  wheel  may  be  Input  to  establish 
the  truck  position,  which  Is  always  assumed  to  lie  parallel  to  the  left 
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and  right  edges  of  the  deck.  The  program  uses  the  stringer  method  with 
each  loaded  element  to  solve  for  the  equivalent  nodal  forces. 

A  set  command  and  regeneration  programs  were  also  written  to 
expedite  the  study.  These  programs  were  used  to  operate  on  the  first 
GTSTRUDL  program  by  rewriting  the  joint  load  Input  portion  based  on  the 
output  of  SKEW  LOADER.  In  this  way  the  truck  could  be  repositioned 
automatically  and  sent  back  to  GTSTRUDL  quasi -Interactively.  This 
provided  for  an  orderly  and  systematic  approach  to  data  acquisition. 
This  program  set  was  written  In  command  language  on  a  VAX  11-750 
computer.  The  program  SKEW  LOADER  is  written  In  FORTRAN-77. 

The  rest  of  Appendix  A  gives  a  flowchart,  program  description  and 
Input  guide  for  SKEW  LOADER.  Copies  of  the  program  are  also  included 
along  with  two  examples  and  sample  outputs. 


A.  2  FLOWCHART 
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INITL 


GNESH 


LOCATE 


00  FOR  EACH  WHEEL 


FIND  WHICH  ELEMENT  THE  WHEEL  FALLS  IN 
BY  SEARCHING  FIRST  IN  THE  SPAN  THEN  IN 
THE  TRANSVERSE  DIRECTION 


COMPUTE  THE  FOUR  NODAL  POINT  REACTIONS 
FOR  THE  ELEMENT  BY  USING 


IlfiMflliWU'M  I  iWMWM 


KEEP  RUNNING  TOTAL  OF  NODAL 
REACTIONS  IN  MATRIX  REACT  (I,  J) 


PRINT  OUT  NON-ZERO  NODAL  POINT 
REACTION  WHICH  6IVES  THE 
EQUIVALENT  DECK  LOADING 


PERFORM  STATIC  CHECK 
TOTAL  NODAL  LOAD  +  WT  OF  OFF  WHEELS 
TRUCK  WT 
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A. 3  Program  Description 


A. 3.1  Main 


The  program  SKEW  LOADER  is  structured  as  a  series  of  subroutine 
calls  from  the  main.  The  initial  matrix  sizes  are  dimensioned 
arbitrarily  In  the  opening  statement  allowing  for  a  maximum  mesh  size  of 
30  nodes  each  way  and  a  maximum  truck  size  of  18  wheels.  These  numbers 
are  only  set  for  convenience  so  that  variably  dimensioned  subroutines 
can  be  used  later.  Therefore  the  maximums  above  can  be  Increased  easily 
by  adjusting  the  sizes  in  the  opening  dimension  statement  (no  other 
program  changes  are  required).  The  Input  is  read  from  L0ADER.DAT  and 
the  output  goes  to  OUTPUT. LIS.  Only  the  key  subroutines  will  be 
discussed  below  as  the  functions  of  the  others  are  self  evident.  The 
subroutine  descriptions  below  contain  a  complete  listing  of  the 
variables  used  within.  Variables  which  are  not  listed  are  defined 
within  the  subroutine  in  terms  of  the  ones  described  here. 


A.3.2  Subroutine  INITL 

This  is  where  most  of  the  user  supplied  data  is  read  In,  stored  and 
echo  printed.  The  following  variables  are  used  here: 


VARIABLE 

AXLSP 


AXLWT 


DDYSK 


DESCRIPTION 
Axle  spacing  matrix 
Axle  weight  matrix 
Deck  dimension  along  X-axis 
Deck  dimension  along  Y  (skew) -axis 
Error  switch 

Mesh  print  request  switch 


V-V.\\  .v- 


y.v* 


i  t 


L  'Hi  ii/jHe*.. 


•v*o*va:**^ 


-> V  \  --V/ V-V 


VARIABLE 

DESCRIPTION 

LCN 

Load  case  number 

NA 

Number  of  axles 

NAXLSP 

Number  of  axle  spaces 

NW 

Number  of  wheels 

NXCROS 

Number  of  crossings  (mesh  lines)  X  direction 

NXSTEP 

Number  of  steps  (elements)  X  direction 

PHI 

Angle  of  deck  skew 

TRTYPE 

Truck  type  (SU-4  or  custom) 

WIDTH 

Width  of  truck 

WJX 

X  coordinate  of  user's  Input  wheel 

WJY 

Y  coordinate  of  user's  Input  wheel 

WNUM 

Input  wheel  number 

WT 

Wheel  weight  number 

XINC 

Increment  (element  size)  X-axis 

XSKINC 

Increment  (element  size)  Y  (skew) -axis 

XSHIFT 

% 

Shift  coordinate  switch  X-axis 

YSHIFT 

Shift  coordinate  switch  Y-axis 

A.3.3  Subroutine  GMESH 

The  mesh  Is  generated  here  based  on  the  user  Inputs  specified  In 
IN ITU  The  nodes  are  generated  first  In  the  span  direction  and  then  along 
the  X-axis.  The  coordinates  of  each  node  are  stored  In  column  vectors  XCOR 


VARIABLE 

ERR 

KEY 

NXCROS 

HXSTEP 

NYCROS 

NYSTEP 

PHI 

XCOR 

XINC 

YCOR 

YSKINC 


DESCRIPTION 
Error  switch 

Mesh  print  request  switch 

Number  of  crossings  (mesh  lines)  X  direction 

Number  of  steps  (elements)  X  direction 

Number  of  crossings  (mesh  lines)  Y  direction 

Number  of  steps  (elements)  Y  direction 

Angle  of  deck  skew 

Nodal  X  coordinate  storage  matrix 

Increment  (element  size)  X-axis 

Nodal  Y  coordinate  storage  matrix 

Increment  (element  size)  Y  (skew) -axis 


A. 3. 4  ^Subroutine  LOCATE 

Now  that  the  position  of  one  of  the  truck  wheels  is  known  (user's 
input)  the  rest  can  be  determined  assuming  that  the  truck  is  lying  along 
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the  span  direction  (i.e.t  parallel  to  the  deck's  free  edges).  From 
trigonometry  the  position  of  the  left  front  wheel  can  be  found  from  the 
known  position  of  the  Input  wheel.  Then  a  standard  procedure  can  be 
followed  to  locate  all  wheels  relative  to  this  position  by  geometry.  It 
may  often  occur  that  one  or  more  of  the  truck  wheels  are  lying  outside  the 
boundary  of  the  deck.  In  this  case  a  warning  Is  printed  to  show  the  user 
that  the  total  weight  of  the  truck  Is  not  on  the  deck.  This  must  be  taken 
Into  account  In  the  statics  check  at  the  end  of  the  program. 


The  following  variables  are  used  in  LOCATE: 


VARIABLE 

AXLSP 

DOX 

DDYSK 

ERR 

NA 

NAXLSP 

NU 

NXCROS 

NYCROS 

OFFWT 

PHI 

PSD 

REACT 

USERX 

USERY~~ 

WHEELX 

WHEELY 

WIDTH 

WJX 

WJY 

WT 

W1X 

W1Y 

XSHIFT 

YSHIFT 


This 


DESCRIPTION 

Axle  spacing  matrix 

Deck  dimension  along  X-axis 

Deck  dimension  along  Y  (skew) -axis 

Error  swtich 

Number  of  axis 

Number  of  axle  spaces 

Number  of  wheels 

Number  of  crossing  (mesh  lines)  X  direction 
Number  of  crossing  (mesh  lines)  Y  direction 
Matrix  containing  weights  of  wheels  off  deck 
Angle  of  deck  skew 
Perpendicular  support  distance 
Matrix  of  nodal  reactions 
X  input  (carried) 

Y  input  (carried) 

X  coordinate  of  wheel  under  consideration 

Y  coordinate  of  wheel  under  consideration 
Truck  width 

X  coordinate  of  user's  input  wheel 

Y  coordinate  of  user’s  input  wheel 
Matrix  of  wheel  weights 

X  coordinate  of  wheel  #1 

Y  coordinate  of  wheel  #1 
Shift  coordinate  switch  X-axis 
Shift  coordiante  switch  Y-axis 

3.5  Subroutine  SEARCH 

subroutine  Is  the  heart  of  the  program  and  certainly  the  most 


complex.  The  function  here  Is  to  find  the  location  of  each  wheel  with 
respect  to  Its  surrounding  nodes,  then  compute  the  nodal  point  reactions  by 
the  stringer  method.  For  a  right  angle  deck  the  search  is  a  relatively 


simple  procedure  as  the  nodal  coordinates  can  be  easily  compared  to  the 
coordl antes  of  the  wheel.  The  skewed  deck  however,  provides  for 
complications  In  locating  the  position  of  the  wheel  relative  to  Its 
neighboring  nodes  as  a  series  of  offsets  from  the  edge  of  the  skewed 
element  must  be  computed  and  compared.  The  idea  is  to  first  search  the 
span  direction  for  the  first  node  with  a  Y  coordinate  larger  than  the  Y 
coordinate  of  the  wheel.  Calling  the  Y  coordinate  of  this  node  and  the  one 
immediately  below  it  D1  and  02,  respectively,  the  lower  "track"  D2  can  then 
be  searched.  Now  the  first  node  along  this  track  with  an  X  coordinate 
larger  than  the  X  coordinate  of  the  wheel  is  rioted.  The  X  coordinate  of 
this  node  is  called  CHECK1.  Now  the  X  offset  distance  of  the  wheel  from 
the  element  edge  is  examined  to  see  on  which  side  of  the  line  it  falls  (see 
Figure  A.l).  If  XOFFI  is  greater  than  WX0FF1  the  wheel  is  in  the  element 
on  thee. right.  If  It  is  less  than  WX0FF1  then  it  is  in  the  element  on  the 
left.  If  the  two  distances  are  equal  then  it  is  on  the  element  edge.  A 
similar  procedure  is  used  to  locate  the  neighboring  nodes  above  the  wheel 
(i.e.,  on  upper  track  01).  Provisions  are  made  for  cases  in  which  the 
wheel  falls  directly  on  a  node  or  on  an  element  edge. 

Now  that  the  coordinates  of  the  four  neighboring  nodes  are  known,  the 
distances  to  the  wheel  can  be  easily  computed.  An  immaglnary  grid  of 
stringers  Is  laid  on  the  element  as  follows: 

1.  One  is  positioned  under  the  wheel  and  laid  horizontally  just 
reaching  the  element  edges 

2.  Then  two  more  are  laid  under  the  first  along  the  skewed  edges 
The  reactions  of  the  four  corners  from  the  weight  of  the  wheel  are  then 
computed  by  statics  and  a  running  total  of  the  results  Is  kept  In  a  storage 
matrix  for  the  mesh  called  REACT.  The  following  variables  are  used  in 
subroutine  SEARCH: 


XINC 


VARIABLE 


DESCRIPTION 


I,M 

Node  counters 

K 

Wheel  counter 

LCOUNT 

Element  counter 

NA,NB,NE,NW 

Nodes  above,  below,  east  and  west 

NNE,NNW,NSE,NSW 

Node  northeast,  northwest,  etc. 

NU 

Number  of  wheels 

NXCROS 

Number  of  crossings  (mesh  lines)  X  direction 

NXSTEP 

Number  of  steps  (elements)  X  direction 

NYCROS 

Number  of  crossings  (mesh  lines)  Y  direction 

NYSTEP 

Number  of  steps  (elements)  Y  direction 

PHI 

Angle  of  deck  skew 

RA.RB.RL.  RR 

Reaction  above,  below,  left,  right 

REACT 

Matrix  of  nodal  reactions 

RNNE ,RNNW ,RNSE ,RNSW 

Reactions  node  northeast,  etc. 

WHEELX 

X  coordinate  of  wheel  under  consideration 

WHEELY 

Y  coordinate  of  wheel  under  consideration 

VIT  ^ 

Matrix  of  wheel  weights 

XCOR 

Nodal  coordinate  storage  matrix 

XINC 

Increment  (element  size)  X-axis 

YCOR 

Nodal  coordinate  storage  matrix 

YSKINC 

Increment  (element  size)  Y  (skew) -axis 

A.3.5  Subroutine  RESULT 

The  non-zero  nodal  reactions  obtained  from  subroutine  SEARCH  are 
printed  out  here  and  summed  to  check  statics.  The  weight  of  the  truck  Is 
checked  against  the  sum  of  the  total  load  on  the  deck  and  the  weight  of  the 
"off  wheels".  The  user  can  verify  this  result  to  be  the  original  weight  of 
the  truck.  The  following  variables  are  used  here: 

VARIABLE  DESCRIPTION 

NCOUNT  Node  counter 

NXCROS  Number  of  crossing  (mesh  lines)  X  direction 

NYCROS  Number  of  crossing  (mesh  lines)  Y  direction 


OFFWT 

Total  weight  of  wheels  off  the  deck 

REACT 

Matrix  of  nodal  reactions 

TOTAL 

Running  sum  of  reactions  for  statics  check 

TRKWT 

Back  summed  weight  of  truck 

USERX 

X  Input  (carried) 

USERY 

Y  input  (carried) 

WJX 

X  coordinate  of  user's  Input  wheel 

VI JY 

Y  coordinate  of  user's  input  wheel 

WNUM 

Input  wheel  number  (user's  choice) 

XSHIFT 

Shift  coordinate  switch  X-axis 

YSHIFT 

Shift  coordinate  switch  Y-axis 

A.4  Input  Guide 


A.4.1  Data  deck 

The  data  file  must  be  titled  LOADER.DAT  and  be  accessible  from  the 
main.  -The  following  data  must  be  Input  In  the  order  shown  beginning  In 
column  #1.  The  term  "card"  refers  to  one  line  of  data  In  the  file. 


CARD 

FORMAT 

DESCRIPTION 

• 

1 

Free 

Alpha-numeric  title  of  output 

2 

111 

Enter  truck  type  as  follows 

(a)  1  Is  FOOT  SU-4 

(b)  9  Is  custom* 

3 

1F10.1 

Angle  of  deck  skew  (degrees) 

4 

112 

'Load  case  number  (user's  option) 

5 

112 

Input  wheel  number 

6 

111 

XSHIFT  switch 

(a)  1  Is  on 

(b)  0  Is  off 


Ill 


YSHIFT  switch 


(a)  1  Is  on 

(b)  0  Is  off 

1F7.2  X  Input  for  wheel  given  on  card  #5  above  as 

follows: 

(a)  if  XSHIFT  Is  on,  give  the  X 
distance  of  the  wheel  from  the 
left  skewed  edge  of  the  deck 

(b)  if  the  XSHIFT  is  off,  give  the  X 
coordinate  of  the  wheel 

1F7.2  Y  Input  for  wheel  given  on  card  #5  above  as 

follows: 

(a)  If  YSHIFT  is  on,  give  the  distance 
to  the  wheel  as  measured  along  the 
inclined  span  edge  for  the  base 
support 

(b)  If  YSHIFT  is  off,  give  the  Y 
coordinate  of  that  wheel 

111  Mesh  print  request  switch 

(a)  1  Is  yes 

(b)  0  Is  no 

112  Number  of  elements  in  the  X  direction 

112  Number  of  elements  along  the  Y  (skewed) 

direction 

Element  size  X  direction 
Element  size  Y  (skewed)  direction 


1F6.2 

1F6.2 


♦Note:  If  the  custom  option  Is  chosen,  then  the  following  data  must  be 
entered  on  the  next  data  cards  (i.e.,  after  card  #2): 


(a)  Number  of  axles  [format  111] 

(b)  Truck  width  [format  1F6.2] 

(c)  Spacing  for  each  successive  axle  as  measured  from  the 
previous  axle  (one  per  data  card)  [format  1F6.2] 

(d)  Axle  weights  (one  per  data  card)  [format  1F6.23 

A.4.2  Examples 

The  above  proceure  Is  Illustrated  by  two  examples.  The  first  Is  an 
Input  to  SKEW  LOADER  for  the  SU-4  case,  and  the  second  is  for  the  custom 
option.  Figures  A.2  through  A.S  given  the  Input  data  and  resulting  truck 
position  diagrams  for  the  two  cases.  Comments  are  provided  on  the  sample 
inputs  for  clarity. 

A.5  Extension  of  the  Program 

The  program  can  be  easily  extended  to  Include  a  number  of  standard 
truck  types  beyond  the  SU-4.  The  easiest  way  would  be  to  write  more 
subroutines  In  the  form  of  TRTYPE1  to  Include  data  for  any  trucks.  Then  of 
course  modify  the  switches  In  subroutine  INITL.  The  program  could  also  be 
easily  extended  to  Include  more  than  one  truck  on  the  deck  at  a  time. 
Finally,  a  more  sophisticated  mech  generator  could  be  included  to  allow  for 
elements  of  different  dimensions.  This  feature  would  be  particularly 
useful  In  the  analysis  decks  with  stiffening  girders. 


W  LOADER  EXAMPLE  1  (SU-4) 


UT  CHECK 


•SUCK  TYPE 

FDOT  SU-4 

:^4GLE  OF  SKEW  (DEG) 

40.0 

|bAD  CASE  NUMBER 

1 

ylPUT  WHEEL  NUMBER 

5 

jj  INPUT  (SHIFT  ON) 

;  INPUT  (SHIFT  ON) 

192.00 

240.00 

f§£SH  PRINT  REQUEST 

NO 

I 

1 


1BER  OF  ELEMENTS 
EMENT  SIZE 
CK  DIMENSION 


X-AXIS 


12 

48.  00 
576. 00 


Y-AXIS  (SKEWED! 


8 

48.  00 
384. 00 


#**WARNING****  -WHEEL  *  1  IS  OFF  THE  DECK  TO  THE  NORTH  ****WARNING**** 
***WARNING***«  WHEEL  #  2  IS  OFF  THE  DECK  TO  THE  NORTH  ****WARNING**** 


****************************** 
****************************** 
*  * 
*  RESULTS  OF  LATEST  ANALYSIS  * 

####****#**####**#*####**#**#* 

****************************** 


WHEEL  # 


WEIGHT 

IN  ELEMENT  * 

AFFECTING  NODES 

X  COORD 

Y  COORD 

9.  35 

RT  OF 

31 

ONLY 

43 

44 

5.59 

222.  15 

9.  35 

56 

62  63 

71 

72 

63.81 

271.00 

9. 35 

ABOVE 

29 

ONLY 

42 

37.73 

183.85 

9.  35 

55 

61  62 

70 

71 

95.95 

232.70 

9.  35 

RT  OF 

28 

ONLY 

40 

41 

69.  87 

145.  55 

9. 35 

54 

60  61 

69 

70 

128.09 

194.  40 

MSULTANT 

NODAL  LOADS 

|de« 

LOAD  (KIPS) 

„  * 

0.39 

r-M 

8.96 

Is 

9.  39 

*43 

8.96 

|»4 

0.39 

60 

6.22 

8.36 

'62 

8.36 

■63 

3.23 

|69 

i. 

0.49 

I70 

0.60 

71 

0.60 

72 

0.23 

BTATIC8  CHECK 

TOTAL  LOAD  ON  DECK  «=  96.10  KIPS 
TRUCK  WEIGHT  *»  7 0.00  KIPS 

****NOTE**** 

SINCE  THE  X SHIFT  IS  ON  X  INPUT  OF  192.00  ©IVES  X  COORD  37.73 


****NOTE**** 

SINCE  THE  YSHIFT  IS  ON  Y  INPUT  OF  240.00  ©IVES  Y  COORD  1S3.83 


input  comments 


[>;ew  loader  example  2  (custom) 


.  *  * v  -  *  . 
i.’-y-: 

t  •*. 

i*-\*  v--  <•  <, 


CHECK 


-jmCK  TYPE  CUSTOM 

>  NUMBER  OF  AXLES  3 

I  TRUCK  AXLE  WIDTH  72.00 

Ij  AXLE  SPACE  #  1  168.00 

AXLE  SPACE  *  2  180.00 

*  AXLE  WEIGHT  «  1  8.00 

bAXLE  WEIGHT  *  2  32.00 

AXLE  WEIGHT  ft  3  32.00 

OF  SKEW  (DEG)  30.0 

CASE  NUMBER  1 

WHEEL  NUMBER  3 

i  INPUT  (SHIFT  ON)  72.00 

\  INPUT  (SHIFT  ON)  192.00 

1ESH  PRINT  REQUEST  NO 


v •>  -  - •  *■* <- - 
\v  v.v:v  v 
,-v  • v.v  ,> 


ESH  DATA 


X-AXIS 


Y-AXIS  (SKEWED) 


NUMBER  OF  ELEMENTS 
LEMENT  SIZE 
CK  DIMENSION 


!  V  V* 


*  *  *. 


36.  00 
432.  00 


24.00 

360.00 


♦♦•WARNING****  WHEEL  ft  2  IS  OFF  THE  DECK  TO  THE  NORTH  ****WARNINO**** 


♦♦♦♦♦♦♦♦a********************* 
•a*#***##*##*###*#***#**###**** 
»  • 

*  RESULTS  OF  LATEST  ANALYSIS  * 

*  ♦ 
****************************** 
*****#**#****#**##*♦*#♦#**#*** 


WHEEL  «  WEIGHT  IN  ELEMENT  *  AFFECTING  NODES  X  COORD 


4. 00  ABOVE  30 


ONLY  48 


->108.00 


Y  COORD 


311.77 


t.  ■*. 


WM 


■  «  t  a  *-»  %  *  j  *  .  V.'j’j  'j  '  *  V  j  V  ^  *  .*  *  .*  .•  ,*  «■  /  .'Ml.'  .*  .*  .*  .*  -  .*  *  .*  »>  ■ 


16.00 

v'. J."  V'-.- 

ABOVE  23 

-T™’"  ^ '•7  "• 'V 

ONLY  41 

.  --  --  *.*’  ■O  'v 

-24.  00 

9 

166. 

16.  00 

70 

74  75  90  91 

38.35 

202. 

16.00 

RT  OF  16 

ONLY  33  34 

66.00 

io. 

16.  00 

63 

67  68  83  84 

128.  35 

46. 

ESULTANT  NODAL  LOADS 
3DE  *  LOAD  (KIPS) 

33  8. 00 

34  8.00 

II  16. 00 

18  4.00 

b7  8. 49 

b8  2.56 

74  2. 96 

73  8. 09 

33  3. 80 

34  1.15 

90  1. 33 

91  3. 62 


i  .  -8. 


v--,* 


STATICS  CHECK 

TOTAL  LOAD  ON  DECK  «  68.00  KIPS 
TRUCK  WEIGHT  **  72.00  KIPS 

****NOTE**** 

SINCE  THE  XSHIFT  IS  ON  X  INPUT  OF  72. 00  GIVES  X  COORD  -24.  OO 
•*»*NOTE**»* 

SINCE  THE  Y8HIFT  IS  ON  Y  INPUT  OF  192.00  GIVES  Y  COORD  166.28 


3  3  £ 


ooooonnonooonfiooooooooorjnofjonoooooonnnooDooorioo 


C  ***************************************************************** 

********  **  **  «*«*««**  **  ** 

********  **  **  ********  **  *•* 

**  **  **  **  **  ** 

********  *****  ******  *•  ** 

********  **  **  ******  **  *  ** 

««  ««  **  *#  **  ***  ** 

********  *«  **  ********  **  *•  «*  ** 

********  **  **  ********  **  ** 

**  ********  **  ******  ********  ****** 

**  ********  ****  *******  ********  ******* 

**  **  **  **  **  **  **  **  **  ** 

**  *»  **  ********  **  **  ******  ******* 

**  **  **  ********  **  »*  ******  ****** 

**  «*  **  #*  **  #*  **  **  **  «* 

*******  ********  **  **  *******  ********  **  ** 

*******  ********  **  **  ******  ********  «*  ** 


THIS  PROGRAM  CALCULATES  THE  SET  OF  EQUIVALENT  NODAL  LOADS  FOR 
ANY  TYPE  OF  TRUCK  MOVING  ACROSS  A  SKEWED  (PARALLELOGRAM)  DECK 
OF  ANY  DIMENSIONS  WITH  ANY  MESH  SIZE.  THE  MESH  IS  GENERATED 
INTERNALLY  AND  CAN  BE  OUTPUT  OPTIONALLY.  THE  USER  MUST  INPUT 
THE  FOLLOWING: 

(1)  ' ALPHA '.  GIVE  ON  THE  FIRST  DATA  CARD  THE  ALPHA-NUMERIC 
TITLE  TO  BE  PRINTED  ON  THE  OUTPUT  LISTING.  C77  CHARACTER 
MAXIMUM] 

(2)  'TRTYPE'.  ENTER  HERE  THE  INTEGER  CORRESPONDING  TO  THE 
TRUCK  TYPE 

FOR  WHICH  THE  ANALYSIS  IS  TO  BE  PERFORMED.  THE  OPTIONS 
ARE  LISTED  BELOW  : 

~~  (A)  1  IS  F DOT  TYPE  SU-4 

<B>  9  IS  CUSTOM 


***NOTE***  IF  THE  CUSTOM  OPTION  IS  CHOSEN.  THEN  THE  FOLLOWING 
DATA  MUST  BE  ENTERED  ON  THE  NEXT  DATA  CARDS  : 

(A)  NUMBER  OF  AXLES.  [INTEGER  UP  TO  91 

(B)  TRUCK  WIDTH  (COAXIAL  WHEEL  SPACING).  [REAL] 

(C)  SPACING  FOR  EACH  SUCCESSIVE  AXLE  AS  MEASURED 
FROM  THE  PRECEEDING  AXLE  (ONE  PER  DATA  CARD). 
NOTE  THAT  THE  #  OF  AXLE  SPACES  IS  ONE  LESS  THAN 
THE  NUMBER  OF  AXLES.  [REAL] 

(D)  AXLE  WEIGHTS  (ONE  PER  DATA  CARD  IN  KIPS).  [REAL! 
ALSO  NOTE  THAT  IF  THE  CUSTOM  OPTION  IS  IN  EFFECT  THEN  THE 
NUMBER  OF  DATA  CARDS  WILL  VARY. 


(3)  'PHI'.  THE  ANGLE  OF  SKEW  (DEGREES).  [REAL] 

(4)  'LCN'.  THE  LOAD  CASE  NUMBER  (ONLY  FOR  USER'S  REFERENCE) 
[INTEGER] 

O)  'WNUM '.  HERE  THE  USER  GIVES  AN  INTEGER  INDICATING  WHICH 
ONE  OF  WHEELS  HE  WILL  LOCATE.  THE  PROGRAM  WILL  AUTO¬ 
MATICALLY  LOCATE  THE  REST  (SEE  SKETCH  FOR  RELATIVE 
POSITIONS  AND  NUMBERING  OF  THE  WHEELS).  NOTE  THAT  THIS 
CARD  COMPLETELY  ESTABLISHES  THE  POSITION  OF  THE  TRUCK 
ASSUMING  THAT  IT  IS  TRAVELLING  PARALLEL  TO  THE  CURB 
(SKEWED  EDGE  OF  THE  DECK).  [INTEGER] 

(6)  'XSHIFT'.  THIS  CARD  ALLOWS  THE  USER  THE  OPTION  OF  GIVING 
THE  X  INPUT  AS  EITHER  SHIFTED  OR  NORMAL.  IF  A  SHIFTED  X 
COORDINATE  DESIRED  THEN  SPECIFY  A  '1'  ON  THIS  CARD  IF 
NOT,  GIVE  A  'O'.  [INTEGER] 

(7)  'YSHIFT '.  THIS  CARD  ALLOWS  THE  USER  THE  OPTION  OF  INPUT- 
ING  EITHER  THE  Y  COORDINATE  OR  INCLINED  SPAN  DISTANCE 
ALONG  THE  SKEWED  EDGE.  IF  SKEWED  Y  COORDINATE  IS  DESIRED 
THEN  SPECIFY  A  '1'  ON  THIS  CARD  .  IF  NOT,  GIVE  A  'O'. 
[INTEGER] 

(5)  IF  CARD  #5  IS  'O'  THEN  GIVE  ON  #7  : 


******************************************************************** 


AD-A145  097  EFFECTS  OF  SKEW  ANGLE  ON  SIMPLE  SPAN  BRIDGE  DECKS  UNDER 
SIMULATED  TRUCK  LOADING(U)  FLORIDA  UNIV  GAINESVILLE 
P  M  KUZIO  1984  N66314-72-A-3029 


2/2 


UNCLASSIFIED 


F/G  13/13 


NL 


93 


C  (A)  THE  X  COORDINATE  OF  THE  WHEEL  NUMBER  GIVEN  ON 

C  CARD  *  4.  UREAL: 

C  IF  CARD  #5  IS  '1 '  THEN  GIVE  ON  #7  : 

C  <B>  THE  X  DISTANCE  OF  THE  WHEEL  FROM  THE  LEFT  SKEWED 

C  EDGE  OF  THE  DECK.  [REAL] 

C  (9)  IF  CARD  #6  IS  'O'  THEN  GIVE  ON  #8  : 

C  (A)  THE  Y  COORDINATE  OF  THE  WHEEL  NUMBER  GIVEN  ON 

C  CARD  #  4.  [REAL] 

C  IF  CARD  #6  IS  '1'  THEN  GIVE  ON  #8  : 

C  (B )  THE  DISTANCE  TO  THE  WHEEL  AS  MEASURED  ALONG  THE 

C  INCLINED 

C  _  SPAN  EDGE  OF  THE  DECK.  [REAL] _ _  __ 

C  ***NOTE***  THE  OPTIONS  OF  SKEWED  COORDINATES  ARE  INDEPENDENT 
C  FOR  THE  X  AND  Y  AXES.  IT  IS  RECOMMENDED  THAT  THE  SKEWED  OPTION 
C  BE  USED  WHEN  PLACING  A  WHEEL  ON  ANY  GRID  LINE  OR  MESH  NODE. 

C  ALTHOUGH  IT  IS  OFTEN  CONVENIENT  AT  OTHER  TIMES  AS  WELL. 

C  CIO)  'KEY'.  IF  KEY**0  THE  COORDINATE  MESH  WILL  NOT  BE  PRINTED 
C  OUT.  IF  KEY  =1  IT  WILL.  [INTEGER] 

C  (11)  'NXSTEP '.  THE  NUMBER  OF  ELEMENTS  IN  THE  X  DIRECTION. 

C  [INTEGER] 

C  (12)  'NYSTEP'.  THE  NUMBER  OF  ELEMETS  ALONG  THE  (SKEWED)  Y 

C  DIRECTION.  [INTEGER] 

C  (13)  'XINC'.  THE  X  INCREMENT  (INCHES).  [REAL] 

C  (14)  'YSKINC '.  THE  INCLINED  Y  INCREMENT  (INCHES).  [REAL] 

C***EXAMPLE*** 

C  PLACE  WHEEL  *  3  ON  NODE  (5.4)  OF  A  30  DEG  SKEWED  DECK.  THE 
C  DECK  MEASURES  540"  X  360"  (X  DIMENSION  AND  SKEWED  SPAN.  RESP) 

C  AND  THE  MESH  SIZE  IS  54"  IN  THE  X  DIRECTION  AND  45"  IN  THE  Y 
C  DIRECTION.  THUS  THERE  ARE  10  INCREMENTS  ALONG  THE  X  AXIS  AND  8 
C  INCREMENTS  ALONG  THE  SKEWED  Y  AXIS  (PLEASE  SEE  SKETCH).  SAY 
C  THE  TRUCK  TYPE  IS  FDOT  SU-4.  WE  WOULD  INPUT  AS  FOLLOWS: 

C 

C  CARD  INPUT  COMMENTS 


FDOT  SU-4  EXAMPLE 
1 

30.0 

1 

3 

1 

1 

162.0 

180.0 


10 

8 

54.0 

45.0 


OUTPUT  LISTING  TITLE 
TRUCK  TYPE  IS  FDOT  SU-4 
SKEW  ANGLE  IS  30  DEGREES 
LOAD  CASE  NUMBER  (USER'S 
OPTION) 

INPUT  WHEEL  IS  #3 
THE  XSHIFT  IS  ON 
THE  YSHIFT  IS  ON 
NODE  (5,4)  IS  162"  FROM  THE 
LEFT  EDGE  OF  THE  DECK 
NODE  (5.4)  IS  4  SPACES  AT 
45"  EACH  (INCLINED  DISTANCE) 
FROM  THE  BASE  OF  THE  DECK 
KEY**0  MESH  WILL  NOT  BE 
PRINTED  OUT 

NUMBER  OF  X  STEPS  IS  lO 
NUMBER  OF  Y  STEPS  IS  8 
X  INCREMENT  IS  54" 

Y  (SKEWED)  INCREMENT  IS  45" 


C****************************************************^******#***** 
DIMENSION  XC0R(30,  30),  YC0R(30,30),  REACT(30,  30 ) , 

*  WHEELX ( 18) »  WHEEL Y ( 18) ,  WT( 18) ,  AXLSP (8) »  AXLWT(9) 

INTEGER  XSHIFT, YSHIFT,  ERR,  WNUM,  TRTYPE 

PHAfiAfTPP  At  PUA«77 

OPEN  ( UN  I  T=  1 ,  F I  LE=  '  [ .  P  AUL  ]  LOADER .  DAT  ' ,  STATUS5*  '  OLD  '  ) 

OPEN  (UNIT=6,  FILE='[.  PAUL]  OUTPUT.  LIS',  STATUS** 'NEW') 

CALL  TITLE1 (ALPHA) 

CALL  INITL ( WJX,  WJY,  PHI,  NXSTEP,  NXCROS,  XINC.  NYSTEP,  NYCROS. 

*  YSKINC, LCN,  KEY,  XSHIFT,  YSHIFT,  WNUM,  AXLSP.  NA, NAXLSP.  NW. 


♦  WIDTH,  TRTYPE,  WT,  DDX» DDYSK ) 

CALL  GMESHCPHI,  NXSTEP,  XINQ,  NXCROS,  XCOR,  NYSTEP.  YSKINC. 

*  NYC R OS,  YCOR,  KEY,  ERR) 

CALL  LOCATE(WJX,  WJY,  W1X,  W1Y,  PHI,  WNUM,  ERR,  OFFWT ,  XSHIFT. 

*  YSHIFT,  USERX,  USERY,  AXLSP,  NA,  NAXLSP,  NU>  WIDTH,  WT,  NXCROS, 

•  NYC R OS, REACT,  MHEELX,  WHEEL Y,  DDX,  DDYSK) 

CALL  TITLE2 

CALL  SEARCH  (PHI,  WT,  REACT,  XCOR,  YCOR.  WHEELX,  WHEEL  Y,  NXSTEP. 
A  NYSTEP,  XINC,  YSKINC, NXCROS, NYCROS.  NW) 

CALL  RESULT  (NXCROS,  NYCROS,  REACT,  OFFWT,  XSHIFT,  YSHIFT, 

A  WNUM,  USERX,  USERY.  WJX,  WJY ) 

910  STOP 

END  _ 

SUBROUTINE  INITL(WJX,  WJY,  PHI.  NXSTEP.  NXCROS,  XINC,  NYSTEP, 

•  NYCROS,  YSKINC,  LCN,  KEY,  XSHIFT,  YSHIFT,  WNUM.  AXLSP.  NA. 

€■■»■■■  n w  iiinifwii  int  ■  ■ wmIm « E.^MTNDDjL  mSkSi  * ******************** 

C  SUB  INITL  READS  IN  THE  INPUT  DATA  AND  SETS  VARIOUS  INITIAL 
C  VALUES  BASED  ON  THIS  INFORMATION. 


DIMENSION  AXLSP (8),  WT ( 18),  AXLWTC9) 

INTEGER  XSHIFT, YSHIFT,  WNUM, TRTYPE.  ERR 
WRITE  (6,310) 

FORMAT  (/, 7X,  ' - 

. - /,  7X.  '  INPUT  CHECK  ',  /.  7X,  * - 

READ* Cl 7315)  TRTYPE* 

FORMAT  (111) 

IF  ( TRTYPE. NE. 1)  00  TO  325 
WRITE  (6,320) 

FORMAT  (/,  7X.  'TRUCK  TYPE',  18X.  'FDOT  SU-4'> 

CALL  TRTYP1 (AXLSP,  NA,  NAXLSP. NW,  WIDTH,  WT) 

GO  TO  360 

IF-CTRTYPE. NE. 9)  GO  TO  360 
CALL  CUSTM1 (NA,  NAXLSP.  NW) 

CALL  CUSTM2C  AXLSP,  NA,  NAXLSP,  NW.  WIDTH,  WT,  AXLWT) 
WRITE  (6,330)  NA,  WIDTH 

FORMAT  ( /.  7X.  'TRUCK  TYPE'.  19X.  'CUSTOM',  //»  14X. 
'NUMBER  OF  AXLES'.  BX,  12.  //.  14X,  'TRUCK  AXLE  WIDTH', 
6X.  F6.2) 

DO  340  I**l,  NAXLSP 
WRITE  (6,335)  I.AXLSP(I) 

FORMAT  (/,  14X.  'AXLE  SPACE  #  ',  II.  8X»  F6.  2> 

CONTINUE 
DO  350  I«1,NA 
WRITE  (6,345)  I.AXLWT(I) 

FORMAT  (/.  14X,  'AXLE  WEIGHT  •  ',  II,  7X,  F6. 2) 

CONTINUE 

READ  (1.362)  PHI 

FORMAT  ( 1F10. 1 ) 

WRITE  (6.364)  PHI 

FORMAT  (/.  7X,  'ANGLE  OF  SKEW  (DEO) ',  11 X,  F4.  1 ) 

READ  (1,366)  LCN 
FORMAT  (112) 

WRITE  (6,368)  LCN 

FORMAT  (/,  7X,  'LOAD  CASE  NUMBER  ',  14X.  12) 

READ  (1.370)  WNUM 
FORMAT  (112) 

WRITE  (6.372)  WNUM 

FORMAT  (/,  7X,  'INPUT  WHEEL  NUMBER '.  12X,  12) 

READ  (1,374)  XSHIFT 
FORMAT  (111) 

READ  (1,376)  YSHIFT 
FORMAT  (111) 

READ  (1,378)  WJX 
FORMAT  C1F7.2) 

READ  (1.380)  WJY 
FORMAT  (1F7.2) 

IF  ( XSHIFT. EQ.O)  GO  TO  383 


WRITE  (6.382)  WJX 

FORMAT  (/,  7X.  ' X  INPUT  (SHIFT  0N)M0X,F7.2) 
GO  TO  386 

WRITE  (6.384)  WJX 

FORMAT  (/,  7X»  'X  INPUT  (SHIFT  OFF)  '.  9X.  F7.  2) 
IF  (YSHIFT.EQ.O)  GO  TO  389 

UR  TTF  lijy 

FORMAT  (/»7X.  'Y  INPUT  (SHIFT  ON)  '»  10X,  F7.  2) 
GO  TO  392 

WRITE  (6.390)  WJY 

FORMAT  (/,  7X»  'Y  INPUT  (SHIFT  OFF) ',  9X,  F7.  2) 

READ  (1,393)  KEY 

FORMAT  (111) 

IF  (KEY. Ed. 0)  GO  TO  396 
WRITE  (6,394) 

FORMAT  (/,  7X,  'MESH  PRINT  REQUEST'.  13X.  'YES') 
GO  TO  400 

WRITE  (6,398) 

FORMAT  (/,7X,  'MESH  PRINT  REQUEST',  12X. 'NO' ) 

WRITE  (6,409) 

FORMAT  (///, 7X, ' - 


6X» ~Y-AXIS~( SKEWED) 

/  - u  , - 

.410 
112) 

.412 
112) 

6,41 
/.  7X 
.420 
1F6. : 

,430 
FA- 


/,  7X,  'MESH  DATA',  16X,  'X-*AXIS 


i5i 

YS 

TE 

P* 

6i 

/. 

43 

7X 

4) 

.  '1 

NX 

ST 

EP 

C  SUB  0ME8H  GENERATES  THE  SKEWED  MESH  AND  ALLOWS  FOR  OPTIONAL 
C  PRINTOUT  FOR  KEY-1. 

DIMENSION  XCOR  ( NYCROS,  NXCROS  ) ,  YCOR  ( NYCROS,  NXCROS  ) 
INTEGER  ERR 

PHI-3. 141 99269399*PHI / 1B0. 

X-0.0 
Y-0.0 
A-0.0 
XSTEP-O.  O 
NCOUNT— O 
ERR— O 

DO  990  J-l. NXCROS 
DO  980  1-1,  NYCROS 
NCOUNT-NCOUNT+ 1 
YCOR(  I,  J)-Y 
XCOR  (I,  J)— X 

IF  (KEY. EQ. 0)  GO  TO  570 
WRITE(6,  960)  NCOUNT,  I,  J,  X,  Y 

960  FORMAT  (/,5X.'N0DE#  '.  13,  9X,  '(I,  J)-(  '»  12,  ',  12.  ')  '»  9) 

•  'X— ',  F7. 2,  3X,  'Y— ',F7.  2) 

970  A— A+YSKINC 


** 


980 


990 


940 


X1«A*SXN(PHI) 

Y1-A*C08(PHI ) 

X-XSTEP-X1 

Y-Yl 

CONTINUE 

A-0.0 

Y-0.0 

XSTEP-XSTEP+XINC 

X-XSTEP 

CONTINUE 

RETURN 

END 

SUBROUTINE  LOCATE<WJX,  WJY.  W1X.  W1Y,  PHI.  WNUM,  ERR.  OFFWT, 
XSHIFT.  YSHIFT.  USERX.  USERY.  AXLSP.  NA.  NAXLSP.  NW.  WIDTH.  WT, 
NXCROS. NYC R OS,  REACT.  WHEEL X.  WHEELY.  DDX.  DDYSK) 


SUB  LOCATE  USES  THE  INPUT  WHEEL  COORDINATES  OF  WJX  AND  WJY  TO 
FIND  THE  COORDINATES  OF  THE  LEFT  FRONT  (i.e.  ,  DRIVER'S)  WHEEL. 
THEN  IT  WILL  LOCATE  THE  COORDINATES  OF  THE  REST  BASED  ON  THIS 
POSITION.  IT  ALSO  CHECKS  TO  BE  CERTAIN  THAT  EACH  WHEEL  FALLS 
ON  THE  DECK  ITSELF  AND  PRINTS  A  WARNING  LIST  FOR  THOSE  WHEELS 
THAT  ARE  NOT  WITHIN  THE  DECK  BOUNDARIES.  THE  TOTAL  WEIGHT  OF 
THE  WHEELS  THAT  ARE  OFF  THE  DECK  WILL  BE  CONSIDERED  IN  THE 
_  STATICS  CHECK  AT  THE  END  OF  THE  PROGRAM. 

*  *  *  ***** 

DIMENSION  WT(NW).  AXLSP  (NAXLSP  > .  WHEELX  <  NW  > .  WHEELY  (NW). 
REACT (NYCROS.  NXCROS) 

INTEGER  WNUM.  ERR.  XSHIFT,  YSHIFT 
IF  (WNUM. LE. NW)  GO  TO  602 
WRITE  (1,600) 


600 


602 


607 


910 


942 


944 


946 


948 


990 


992 


FORMAT  (//. 7X.  '***  ERROR  ***  ILLEGAL  WHEEL  NUMBER  INPUT') 
ERR-1 

GO  TO  1090 
USERX-WJX 
USERY-WJY 

IF  (YSHIFT. EQ.O)  GO  TO  607 
WJY«WJY*COS(PHI ) 

IF  (XSHIFT. EQ.O)  GO  TO  910  _ ... 

WJX-WJX-(WJY*SIN<PHI ) /SIN( 1.  970796327-PHI ) > 

Z 1 -WIDTH*COS ( PHI ) 

Z2«WIDTH*8IN<PHI ) 

IF  (WNUM.NE.  1)  GO  TO  942 

W1X-WJX 

W1Y-WJY 

GO  TO  980 

IF  (WNUM.  NE.  2)  GO  TO  944 
W1X-WJX-Z1 
W1Y-WJY-Z2 
GO  TO  980 
A-AXL8P(1) 

IF  (WNUM.  NE.  3)  00  TO  946 
W1X-WJX-(A)*SXN(PHX) 

WlY-WJY+( A)*COS(PHI ) 

00  TO  980 

IF  (WNUM. NE. 4)  GO  TO  948 
WlX"WJX-( ( A)*SIN(PHI )+Zl ) 

W1 Y-WJY+ ( (A) *COS ( PHI ) -Z2 ) 

GO  TO  980 
B»AXL8P(2) 

IF  (WNUM.NE. 9)  GO  TO  990 
W1X«WJX-(A+B)»SIN(PHI) 

W1Y«WJY+(A+B)*C0S(PHI ) 

GO  TO  980 

IF  (WNUM.NE. 6)  00  TO  992 
WlX-WJX-( (A+B)*SIN(PHI )+Zl ) 

W1Y-WJY*( ( A+B ) *COS ( PHI ) -Z2 ) 

GO  TO  980 
C-AXL8PO) 

IF  (WNUM.NE.  7)  GO  TO  994 


******** 


W1X-WJX-(A+B+C)*SIN(PHI) 

W1 Y-WJY+  < A+B+C ) *COS ( PHI ) 

GO  TO  980 

IP  (WNUM.NE.  8)  GO  TO  956 
W1 X»WUX-( ( A+B+C ) *S INC PHI ) +Z 1 ) 

WlY«WJY+(  (A+B+C )*CQS( PHI )-Z2) 

GO  TO  980 
D— AXLSP (4) 

IF  (WNUM.NE.9)  GO  TO  958 
W1X-WJX-(A+B+C+D)*SIN(PHI > 

Ml Y— WJY+  (  A+B+C+D  >  *COS  <  PHI > 

GO  TO  980 

IF  (WNUM.NE. 10)  GO  TO  960 
W1X-WJX~( (A+B+C+D )*SIN( PHI >+Zl ) 

WlY-WJY+(  (A+B+C+D )*COS( PHI  >-Z2) 

00  TO  980 
E-AXLSP(5> 

IF  (WNUM.NE. 11)  GO  TO  962 
W1 X«W JX- ( A+B+C+D+E ) *S I N  <  PH I ) 

W 1 Y=WJY+ ( A+B+C+D+E  >  *C08 ( PHI > 

GO  TO  980 

IF  (WNUM.NE.  12)  GO  TO  964 
W1X-WUX~( ( A+B+C+D+E) *8IN (PHI )+Zl ) 

WlY«WJY+( ( A+B+C+D+E )*COS( PH I )-Z2> 

GO  TO  980 
F-AXLSP(6> 

IF  (WNUM.NE. 13)  GO  TO  966 
W1X«WJX-(A+B+C+D+E+F)*SIN(PHI ) 

W1 Y»WJY+(A+B+C+D+E+F)*COS(PHI ) 

GO  TO  980 

IF  (WNUM.NE. 14)  GO  TO  968 
WlX«WUX-( (A+B+C+D+E+F)*SIN(PHI )+Zl ) 

W 1 Y-WJY+ < ( A+B+C+D+E+F ) *CQS ( PHI ) -Z2  > 

CP  TO  980 

Ok  A  VI  apt  7» 

IF  (WNUM.NE. 15)  GO  TO  970 

W1 X— WUX—  <  A+B+C+D+E+F+G ) *S I N ( PHI ) 

W1Y«WJY+(A+B+C+D+E+F+C)*C0S(PHI ) 

00  TO  980 

IF  (WNUM.NE. 16)  GO  TO  972 
W1 X-WJX- ( ( A+B+C+D+E+F+G ) *SIN(PHI ) +Z 1 ) 

WlY«WJY+( (A+B+C+D+E+F+C)*COS(PHI )-Z2) 

GO  TO  980 
H«AXLSP(8) 

IF  (WNUM.NE.  17)  CO  TO  974 
W 1 X— WUX—  <  A+B+C+D+E+F +0+H ) *8 IN  <  PHI > 

W1 Y-WJY+ ( A+B+C+D+E+F+G+H ) *COS ( PHI ) 

GO  TO  980 

IF  (WNUM.NE. 18)  GO  TO  976 
W1 X-WJX- ( ( A+B+C+D+E+F+G+H ) *8 I N ( PHI > +Z 1 ) 

W 1  Y— WJY+  (  <  A+B+C+D+E+F+G+H  )  *COS  (  PHI  )-Z2> 

00  TO  980 
WRITE  (6.978) 

FORMAT  (//.  7Xi  '***#ERROR***+  ILLEGAL  WHEEL  NUMBER  INPUT 
'  ON  DATA  CARD  •  10') 

ERR-1 

GO  TO  1090 
1-0 

WHEELX ( 1 )— W1X 
WHEELX(2)«W1X+Z1 
WHEELY( 1 )»W1Y 
WHEELY(2)— W1Y+Z2 
I-I+2 

IF  (I.EQ.NW)  00  TO  1010 

a— axlsp ( 1 ) 

WHEELX (3)*W1 X+( A )+SIN (PHI ) 

WHEELX (4) -WHEELX (3)+Zl 
WHEELYO)— W1Y-A*C08(PHI ) 

WHEEL Y ( 4 ) -WHEELY ( 3 ) +Z2 


1010 


1013 


1019 


1*1+2 

IF  (X.EQ.NW)  GO  TO  1010 
B*AXLSP(2) 

WHEELX(9)“W1X+IA+B)*SIN<PHI) 

WHEELX (6)“WHEELX <  5)+Zl 
WHEELY ( 9 ) -W 1 Y- ( A+B ) *C08 ( PH I ) 
WHEELY(6)-WHEELY(9)+Z2 
I- 1+2 

XF  (X.EQ.NW)  GO  TO  1010 
OAXLSP<3> 

WHEELX(7)-W1X+<A+B+C)*SIN(PHI) 

WHEELX C 8 ) -WHEELX < 7 > +Z 1 
WHEELY ( 7 ) “W 1 Y- < A+B +C ) *C08 < PH I ) 

WHEELY < 8 ) “WHEEL Y ( 7 > +Z2 
I-I+2 

IF  (I.EQ.NW)  GO  TO  1010 
D*AXLSP(4) 

WHEELX ( 9 )-WlX+( A+B+C+D) *SIN( PHI ) 

WHEELX ( 1 O > -WHEEL X C 9 ) +Z 1 
WHEELY (9 )*W1Y-<A+B+C+D)*C0S (PHI > 

WHEELY < 10)=WHEELY(9)+Z2 
I-I+2 

IF  (I.EQ.NW)  GO  TO  1010 
E-AXLSPO) 

WHEELX ( 1 1 ) “W1 X+ ( A+B+C+D+E ) *S IN ( PHI ) 

WHEELX < 12 ) “WHEELX ( 1 1 ) +Z 1 

WHEELY ( 1 1 )»W1 Y-( A+B+C+D+E )«COS< PHI ) 

WHEELY ( 12 ) “WHEELY ( 1 1 >+Z2 
1*1+2 

IF  (I.EQ.NW)  GO  TO  1010 
F-AXLSP(A) 

WHEELX ( 13 ) *W1 X+ ( A+B+C+D+E+F ) *SXN ( PHI ) 

WHEELX ( 14 ) “WHEELX  < 1 3 ) +Z 1 

WHEELY ( 13)=WlY-( A+B+C+D+E+F ) *COS ( PH  X  > 

WHEELY ( 14 ) “WHEELY ( 13) +Z2 
1*1+2 

IF  (X.EQ.NW)  GO  TO  1010 
C*AXLSP(7> 

WHEELX ( 1 9 ) =W1 X+  < A+B+C+D+E+F+G ) *SIN ( PHI > 

WHEELX ( 16 ) “WHEELX ( 1 9 ) +Z 1 

WHEELY ( 1 5 ) “W1 Y- < A+B+C+D+E+F+G ) *COS ( PHI ) 

WHEELY (16) “WHEELY (19) +Z2 
1*1+2 

IF  (I.EQ.NW)  GO  TO  1010 
H*AXL8P(8) 

WHEELX ( 1 7 ) “W1 X+ ( A+B+C+D+E+F+G+H ) *8IN(PHI ) 
WHEELX ( 1 8 ) -WHEELX < 1 7 ) +Z 1 

WHEELY ( 1 7 ) *W1 Y- ( A+B+C+D+E+F+G+H ) *COS (PHI ) 
WHEELY ( 18) “WHEELY ( 17)+Z2 
1-1+2 

IF  (I.EQ.NW)  GO  TO  1010 
WRITE  (6.990) 

FORMAT  (//. 7X.  'ILLEGAL  NUMBER  OF  WHEELS  INPUT' 
GO  TO  1090 
OFFWT-O. O 
DO  1080  1*1.8 
A*WHEELY(I) 

IF  (A.GE.0.0)  GO  TO  1019 
OFFWT  *OFFWT +WT4  X ) 

WT( I )— 0. 0 
WRITE  (6.  1013)  I 

FORMAT  (///,  7X»  '****WARNING****  WHEEL  *  ',  12, 

'  XS  OFF  THE  DECK  TO  THE  SOUTH') 

00  TO  1080 
P8D—DDY8K+C0S ( PH  X ) 

XF  (A.LE.P8D)  00  TO  1029 
OFFWT “OFFWT +WT ( I ) 

WT  ( I )  “0.  O 
WRITE  (6, 1020)  I 


1020  FORMAT  (///,  7X,  '♦^♦JWARNINC****  WHEEL  •  '.12. 

*  '  18  OFF  THE  DECK  TO  THE  NORTH') 

GO  TO  1080 

1023  OFFSET-A*SIN(PHI > /SIN( 1. 370796327-PHI ) 

B-WHEELX<I> 

Bl-DDX-OFFSET  . 

IF  (B.LE.Bl)  00  TO  1033 
OFFWT-OFFWT+WT ( I > 

WT(I>-0.0 

1030  FORMAT  i///???? '♦♦♦♦WARNING****  WHEEL  *  '. 12. 

*  'IS  OFF  THE  DECK  TO  THE  EAST') 

GO  TO  1080 

1033  B  1-0.0— OFFSET 

IF  (B.GE.B1)  GO  TO  1080 
OFFWT — OFFWT +WT ( I ) 

WT(I>«0.0 
URITE  (A,  1040)  I 

1040  FORMAT  (///,  7X»  '****WARNINO****  WHEEL  #  '.12. 

*  '  IS  OFF  THE  DECK  TO  THE  WEST') 

1080  CONTINUE  _ 

640  DO  630  I— 1»  NYCROS 
DO  630  J-l.  NXCROS 
REACT (I,  J)*=0.  O 
650  CONTINUE 
1090  RETURN 

run  _ 

SUBROUTINE  SEARCH < PHI.  WT.  REACT,  XCOR.  YCOR.  WHEELX. 

*  WHEELY,  NXSTEP.  NYSTEP.  XINC.  YSKINC,  NXCROS,  NYCROS,  NW) 
C********“**i*********S***************«^****^*******?***** 
C  SUB  SEARCH  IS  USED  TO  FIND  THE  LOCATION  OF  EACH  WHEEL  WITH 

C  RESPECT  TO  ITS  SURROUNDING  NODES.  DO  THIS  5XrrFi5iT^S^SnT»jlTCo 
C  VERTICALLY  TO  GET  DISTANCES  D1  &  D2  WHICH  ARE  THE  COORDINATES 
C  OF  THE  "NODES  IMMEDIATELY  ABOVE  8t  BELOW  THE  WWEQ.,  RESP.  THEN 
C  SEARCH  HORIZONTALLY  ALONG  THE  ^COORDINATE  0FD2  TO  OET  THE 
C  DISTANCES  D3  AND  D4  WHICH  ARE  THE  PROPER  X  COORDINATES  OF  THE 
C  NODES  ALONG  THIS  LOWER  'TRACK'.  THEN  USE  A  SIMILAR  PROCEDURE 

C  FOR  D3  AND  D6  ALONG  THE  UPPER  'TRACK'.  _ 

C  ******************************************************** *  *  *  •*** 
DIMENSION  WT(NW),  REACT (NYCROS,  NXCROS). 

*  XCOR (NYCROS,  NXCROS ) .  YCOR  (NYCROS,  NXCROS) ,  WHEELX  ( NW> , 

*  WHEELY(NW) 

DO  830  K-l.NW 

WTCHK— WT(K )  _ _ _ 

IF  (WTCHK. EQ. 0.0)  GO  TO  830 
1-1 

LCOUNT— O 

703  C— YCOR(  I.  1 ) 

D-WHEELY(K) 

CALL  RNDOFF(C.D)  _ 

IF  (C.EQ.D)  GO  TO  710 
IF  (C.GT.D)  GO  TO  720 
1-1*1 

LCOUNT -LCOUNT + 1 
GO  TO  703 
710  11-1  „ 

Dl-YCORdl.  1) 

00  TO  721 


720  I l-I 

12— I 1—1 

NNE-I1 

NSE-I2 

Dl-YCORdl.  1) 
D2-YC0R( 12,  1) 

GO  TO  729 

721  N«1 

LCOUNT -LC0UNT-NY3TEP 
F— WHEELX(K) 


********** 


C— XCOR ( II*  N>  _  % 

CHECKS* AB8 < WHEELX ( K ) - < N- 1 > *X I NC ) 

ChIcK4*AB8(WHEELY(K>*8IN<PHI)/SINC1.  970796327-PHI)  > 
CALL  RNDOFF < CHECK3,  CHECK4) 

CALL  RNDOFF <G*  F)  _ 

IF  (CHECK3.EQ.CHECK4)  00  TO  723 
IF  (O.OT.F)  CO  TO  729 


LcSuNT-LCOUNT+NYSTEP 


'ONLY 


NE-NE+NYCROS 
CO  TO  722 
ISM*N 

D9-XC0R(U*  19) 

W-WT(K) 

A-WHEELX(K) 

BAMEELY(K) 

WRITE  (6.724)  K,  W*  LCOUNT.  NE*  AjB  _ w  , 

FORMAT  ( /*  9X,  1 1#  6X#  F9.  2*  4X*  'ABOVE  '»  12*  8X»  'ONLY 
12*  7X*  F7.  2*  2X*  F7.  2) 

RS-WT(K) 

REACT (II*  I9)»REACT(11.  I9)+RS 
00  TO  890 
I9=N 
I6-N-1 

D9-XC0R(U* 19) 

D6»XC0R(I1.  16) 

NW— NE— NYCROS 
W-WT(K) 

A-WHEELX(K) 

B-WHEELY(K)  .  _ 

WRITE  (6*  726)  K,  W*  LCOUNT*  NW.  NE,  A.  B 
FORMAT  (/,  9X*  I1*6X*F9.2*  4X*  'ABOVE  '*  12*  8X.  'ONLY 
12*  IX*  12.  4X*  F7.  2.  2X.  F7.  2) 

RL- ( D9-A ) *WT < K ) / ( D9-D6 ) 

RR— WT(K)— RL 

REACT (II* I9)“REACT( II* I9)+RR 
REACT (II. 16)— REACT ( II* I6)+RL 
00  TO  890 
L-l 

LCOUNT-LCOUNT-NYSTEP 

E-XC0R(I2*L) 

i  hiti  x  r »r  i 

CHECK4— ABS  ( WHEEL Y ( K ) *8 IN ( PHI )  /S IN ( 1  •  970796327-PHI )  > 
CHECK9-AB8 ( WHEELX (K)-( L-l )*X INC) 

CALL  RNDOFF ( CHECK9.  CHECK4 > 

CALL  RNDOFF(E.F) 

IF  (CHECK9. EQ.  CHECK4)  00  TO  731 

IF  (E.OT.F)  00  TO  740 

L-L+l 

LCOUNT-LCOUNT+NYSTEP 
N6E-N8E+NYCR08 
00  TO  730 
I3-L 
19-13 

D'3— XCOR  (12*  13) 

D9—XC0R (II* 13) 

NB-N8E 

NA-NB+1 

W-WT(K) 

A-WHEELX(K) 

B -WHEEL Y ( K ) 

WRITE  (6,732)  K,  W*  LCOUNT,  NB.  NA,  A,  B 


FORMAT  (/,  9X,  1 1 ,  6X,  F9.  2,  4X,  'R 
12,  IX,  12,  4X,  F7.  2,  2X,  F7.  2) 

RA— WT(K)*(B-D2)/(D1— D2) 
RB-WT(K)-RA 

REACT ( 12,  13) —REACT (12,  I3)+RB 
REACT (II, I9)-REACT(I1,  I9)+RA 
00  TO  890 


RT  OF  12,  8X*  'ONLY 


740 


730 

760 


770 


X0FF1«(D-D2)*SIN(PHI >/SIN< 1. 570796327-PHI > 
CHECK1"XC0R  <  12.  L) 

WX0FF1 "CHECK 1-F 

CALL  RNDOFF (WXOFFl.  X0FF1 > 

IF  (WXOFFl.  GT.X0FF1)  00  TO  730 

13- L+1 
I4"L 

LCOUNT"LCOUNT +N YSTEP 
NSE-NSE+NYCROS 
00  TO  760 
I3"L 

14- L-1 

D3"XC0R( 12. 13) 

D4»XC0R(I2> 14) 

NSW^MSE— NYCROS 

N-l 

G"XC0R(I1.N) 

CALL  RNDOFF(G.  F) 

IF  (G.CT.F)  GO  TO  780 


780 


790 


793 


NNE-NNE+NYCROS 
00  TO  770 

XOFF2«(Dl-D)*SIN(PHI > /8INC 1 . 570796327-PHI  > 

CHECK2=XC0R ( 1 1,  N) 

WX0FF2»F- C  CHECK2-X I NC ) 

CALL  RNDOFF  <  WX0FF2. X0FF2) 

IF  (WX0FF2.  OT.  X0FF2)  00  TO  790 
I3-N-1 
I6-N-2 

NNE-NNE-NYCROS 
LCOUNT "LCOUNT— N YSTEP 
OO  TO  795 
I3-N 
I6*N— 1 

D3"XC0R(I1. 13) 

06"XC0R (Il« 16) 

NNW-NNE— NYCROS 
W-WT(K) 

A"UHEELX(K> 

B-WHEELY(K) 

WRITE  (6.797)  K.  W.  LCOUNT.  NSW.  NNW.  NSE.  NNE.  A.  B 
FORMAT  (/.  9X.  II.  6X.  F3.  2.  7X.  12.  10X.  12.  IX.  12.  IX.  12. 

•  IX.  I2,4X»F7.2.2X,F7.2) 

C  COMPUTE  NODAL  POINT  REACTIONS  BY  PUTTING  AN  IMAGINARY  STRINGER 
HORIZONTALLY  ACROSS  THE  ELEMENT  AND  FINDING  THE  EQUIVALENT 
STATIC  LOADS  ON  THE  SKEWED  EDGES.  THEN  PUT  TWO  IMAGINARY 
STRINGERS  ALONG  THE  SKEWED  EDGES  AND  SOLVE  FOR  THE  NODAL  POINT 

LOADS.  AS  THE  REACTIONS  FROM  EACH  WHEEL  ARE  FOUND.  KEEP  A _ 

RUNNING  TOTAL  OF  EACH  WHEEL'S  CONTRIBUTION  IN  A  STORAGE  MATRIX 
REACT(I.  J). 


797 


800 


850 


D7*WHEELX(K>— D4 
DB*  ( WHEEL  Y  (  K  >  -D2 )  *S  I N  ( PH  I )  /SINC 1.  370796327-PHI ) 
RL-WT  C  K ) * ( <  D3-D4 ) - ( D7+D8 ) > / ( D3-D4 ) 

RR»WT(K)-RL 

RNNW-RL*  <  WHEEL Y  <  K ) -D2 ) / <  D 1 -D2 ) 

RN8W"RL* ( D 1 -WHEEL Y ( K ))/( D 1 -D2 ) 

RNNE"RR* ( WHEEL Y ( K ) -D2 )/< D 1 -D2 ) 

RN8E-RR*  <  D 1 -WHEEL Y (  K )  >  /  ( D 1 -D2 ) 

REACT (II. I6)"REACT(I1.  I6)+RNNW 
REACT (12,  14 ) "REACT ( 12.  I4)+RNSW 
REACT (II. I3)"REACT( II. I5)+RNNE 
I3)-REACT(I2,  I3)+RNSE 


913 


REACT  (12. 

CONTINUE 
RETURN 
END 

SUBROUTINE  RNDOFF ( X.  Y) 


€»»»»»»»»»»»»»«»»»»»»•»»»»»»»»»»»»*»»»*»»»*»♦♦»»»»»♦»♦♦♦»»»»»»»«•♦♦ 


******* 


C  SUBROUTINE  RNDOFF  ROUNDS  OFF  THE  VALUES  ASSIGNED  TO  X  AND  Y  * 
~  TO  THREE  DECIMAL  PLACES.  .  * 

X«1000.0#X  #*  #  *************  *** 

Y-IOOO. 0*Y 

X*ANINT(X)/1000.0 

Y«ANINT(Y)/1000.0 

RETURN 

END 

SUBROUTINE  RESULT (NXCROS*  NYCROS,  REACT.  OFFWT.  XSHIFT, 

*  YSHIFT.  WNUM,  USERX.  USERY.  WJX.  WJY  > 

C »♦»»»»»»»»»»«»»»»♦»»»»»»»»»»»»»»»»»»»»»»  *****  WHHHHt »»»»»»»»»»«»» 

C  SUB  RESULT  OUTPUTS  THE  RESULTANT  (NON— ZERO)  NODAL  LOADS  AND  » 

C  CHECKS  BY  STATICS  TO  ASSURE  THAT  THE  SUM  OF  THE  WEIGHTS  OF  * 

C  THE  'OFF-WHEELS'  ADDS  UP  TO  THE  TRUCK  WEIGHT.  * 

C<hhhhmmhmmhhmmhhmhhmmmmhmmmhmhhmhhhhhhhhhhhhhhmhhhhhmhhhmhhmhhmm> 
DIMENSION  REACT (NYCROS.  NXCROS) 

INTEGER  XSHIFT, YSHIFT, WNUM 
TOTAL-O.  O 
NCOUNT-1 
WRITE  (6,855) 

855  FORMAT  (  '1  ',////,  7X,  'RESULTANT  NODAL  LOADS', //,  7X,  'NODE  #'. 

*  5X»  'LOAD  (KIPS)  '.  /) 

DO  880  J*l,  NXCROS 
DO  870  1*1,  NYCROS 
R*REACT(I.J) 

IF  (R.EQ.O.O)  00  TO  860 
TOTAL*TOT  AL+R 
WRITE  (6,857)  NCQUNT , R 
857  FORMAT  (/.  8X,  12.  9X.  F7.  2) 

860  NCOUNT  *NCOUNT  + 1 

870  CONTINUE 
880  CONTINUE 

TRKWT “TOT AL+OFF WT 
WRITE  (6,890)  TOTAL,  TRKWT 
890  FORMAT  (////,  7X.  'STATICS  CHECK', //,  7X, 

*  'TOTAL  LOAD  ON  DECK  “  ',F5.2,  '  KIPS'.  //.  7X. 

*  'TRUCK  WEIGHT  =  '.F5.2,  '  KIPS') 

IF  (XSHIFT. EQ.O)  GO  TO  900 
WRITE  (6,895)  USERX.  WJX 

895  FORMAT  (//,  7X,  ' ****NOTE**** ' ,  //,  7X,  'SINCE  THE  XSHIFT'. 

*  '  IS  ON  X  INPUT  OF ',  F7. 2,  '  GIVES  X  COORD  '.F7.2) 

900  IF  (YSHIFT.  EQ.O)  GO  TO  910 

WRITE  (6.905)  USERY.  WJY 

905  FORMAT  (//,  7X.  '****NOTE**** ',//,  7X, 'SINCE  THE  YSHIFT'. 

*  '  IS  ON  Y  INPUT  OF',  F7. 2,  '  GIVES  Y  COORD  '.F7.2) 

910  RETURN 

END 

SUBROUTINE  TRTYP1 (AXLSP,  NA,  NAXLSP,  NW,  WIDTH, WT) 
C***************************************************************** 
C  SUB  TRTYP1  HOLDS  THE  REQUIRED  INPUT  DATA  FOR  THE  FDOT  SU-4  * 
C  TYPE  TRUCK  _  _ * 

DIMENSION  AXLSP (3),  WT(8) 

NA*4 

NAXLSP-3 

NW*8 

WIDTH*76.  0 
AXLSP ( 1 )*1 10. O 
AXLSP (2) *50.0 
AXLSP (3) *50.0 
WT( 1 )*6. 95 
WT(2)«6.  95 
DO  950  1*3,8 
WT ( I ) *9.  35 
950  CONTINUE 
RETURN 
END 

SUBROUTINE  CUSTM1 (NA,  NAXLSP, NW) 


m 


103 


C  »#*#*#*####*###■**■#••»•-*  *****##**#**#■#  #**#***«--*******- «■•#•■«•■»■*••»(-•»■#■»■  •*•**■** 
C  SUB  CUSTM1  READS  IN  THE  NUMBER  OF  AXLES  AND  FINDS  THE  NUMBER  * 
C  OF  AXLE  SPACES  AND  THE  NUMBER  OF  WHEELS.  THE  INFORMATION  IS  * 
C  REQUIRED  FOR  THE  VARIABLE  DIMENSIONS  IN  SUB  CUSTM2.  * 

C****#*********************#********#*********#*******#***#*-***** 
READ  (1.705)  NA 
705  FORMAT  (111) 

NAXLSP»NA-1 

NW*2«NA 

RETURN 

END 

SUBROUTINE  CUSTM2(  AXLSP,  NA.  NAXLSP,  NW»  WIDTH,  WT,  AXLWT) 
C*******************#********#***********#**********'**  «-#******#** 
C  SUB  CUSTM2  READS  IN  THE  REST  OF  THE  REQUIRED  DATA  FOR  ANY  * 

C  TYPE  OF  TRUCK  * 

C*****#***»*********»********#***i»*******i»*«*********tt*»'»'IHHHHHHt'* 

DIMENSION  AXLSP  (NAXLSP),  WT(NUI),  AXLWT  (NA) 

READ  (1,740)  WIDTH 
740  FORMAT  (1F6.2) 

DO  770  1=1, NAXLSP 
READ  (1,750)  AXLSP  ( I  ) 

750  FORMAT  (1F6.2) 

770  CONTINUE 

DO  780  1=1. NA 

READ  (1,775)  AXLWT(I) 

775  FORMAT  (1F6.2) 

780  CONTINUE 

DO  790  K=1.NA 

I»2*K-1 

J=2*K 

WT(I)=AXLWT(K)/2 

WT(J)=WT(I) 

790  CONTINUE 
RETURN 
END 

SUBROUTINE  TITLE1 (ALPHA) 

C#*************#***************«**««************4MHHI'****'iHHHHHHHHt- 

C  SUB  TITLE1  READS  IN  AND  ECHO  PRINTS  THE  TITLE  OF  THE  OUTPUT  * 
C  LISTING  « 

C**************************************************************** 
CHARACTER  ALPHA*77 
READ  (1,200)  ALPHA 
200  FORMAT  (A77) 

WRITE  (6,210)  ALPHA 
210  FORMAT  (///.7X.A77) 

RETURN 

END 

SUBROUTINE  TITLE2 
WRITE  (6.700) 

700  FORMAT  (///.  20X.  '***##******#****##***********«•'./,  20X» 

'#****#*#******#*#****##*****#*',  /,  20X»  '*  ',  28X.  '*'»  /,  20X» 

'*  RESULTS  OF  LATEST  ANALYSIS  *',/.20X.  '*',28X»  '*',/,20X, 
'****************#*************',  /,  20X. 
'**#********»**•*****»*•*******',  ///,  7X»  'WHEEL  #', 

2X,  'WEIGHT',  2X.  'IN  ELEMENT  #',2X.  'AFFECTING  NODES',  2X. 


9 

* 

* 

* 


'X  COORD'. 2X,  'Y  COORD',/) 
RETURN 
END 


TT^T*T^TT*! 

•  ^  f 


l  .  -• 


L.  . 


U*-  - 


t _ 1L 


C****************tt*tt********#******#****^****M*«"^**-fr*'iHHHHt4HMMt 

C  CREATE.  FOR  * 

C*******»##***#*******************#*****#«*-*****-»**#«"*-»****'iHHHHMt 
CHARACTER  BUFFER*9,  BETA*79,  ALPHA*44»  GAMMA*15 
OPEN  (UNIT*! ,  FILE* ' C .  PAUL! OUTPUT.  LIS'.  STATUS* ' OLD ' ) 

OPEN  (UNIT*2,  FILE- '  t.  PAUL3PS20100.  DAT'.  STATUS*  ' OLD '  > 

OPEN  (UNIT=3.  FILE»'C.  PAUL3PS20100.DAT',  STATUS* 'NEW') 

DO  10  1*1. 1000 
READ  (1,5)  BUFFER 
5  FORMAT  (3X»  A9) 

IF  (BUFFER.  EG.  'RESULTANT')  GO  TO  15 
10  CONTINUE 

15  READ  (1.20)  NJ,  FJ 

20  FORMAT  (////.  3X.  13,  11X,  F5.  2) 

FJ*-FU 

READ  (2, ' ( A79 ) ' )  BETA 
WRITE  ( 3,  ' ( A79 ) ' >  BETA 

READ  (2.  '  (  A44.  1F5.  1,  A15)  '  )  ALPHA,  POSIT,  GAMMA 
P0SIT-P0SIT+12.0 

WRITE  (3.  '  (  A44,  F5.  1.  A15)  '  )  ALPHA.  POSIT.  GAMMA 
DO  30  1*1,23 
READ  (2,  ' ( A79) ' )  BETA 
WRITE  (3. ' ( A79) ' )  BETA 
30  CONTINUE 

WRITE  (3,40)  NJ.  FJ 
DO  50  1*1,25 
READ  (1.35)  NJ,  FJ 
35  .  FORMAT  ( /.  3X,  13,  1 IX.  F5.  2) 

FJ*— FJ 

IF-fNJ.EQ.  O)  GO  TO  70 
WRITE  (3,40)  NJ,  FJ 

40  FORMAT  ( 13.  IX.  'FORCE  Z  '.F5.2) 

50  CONTINUE 

70  WRITE  (3,80) 

80  FORMAT  ( 'STIFFNESS  ANALYSIS ', /, 

*  'CALCULATE  AVERAGE  PRINCIPAL  BENDING  RESULTANTS',  /, 

*  'FINISH') 

STOP 

END 


CTSTR.COM 


*  GTSTRUDL  C. PAUL3PS20100. DATj 
«  ©E. PAUL3RUN 

*  GTSTRUDL  C . PAUL3PS20100.  DAT; 

*  ©E. PAUL3RUN 

*  GTSTRUDL  E. PAUL3PS20100. DATj 

*  ©E. PAUL3RUN 

*  GTSTRUDL  C. PAUL3PS20100. DATj 

*  ©E. PAUL] RUN 

*  GTSTRUDL  C.PAUL3PS20100.DAT! 

*  @E. PAUL DR UN 
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FIGURE  C.5  CONTOUR:  Z  DISPLACEMENT  FOR  CENTRAL  10  KIP  LOAD,  SKEW 


FIGURE  C.8  CONTOUR:  Z  DISPLACEMENT  FOR  PEAK  SU-4 
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